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Introduction 

One  of  the  primary  objectives  of  this  report  is  to  provide  an  overview  of  the  sources  of  noise 
in  engine  test  cells.  An  engine  test  cell  is  either  a  closed  or  a  partially  closed  system.  The  solid 
surfaces  that  form  the  cell  reflect  back  any  incident  acoustic  waves.  Thus  a  test  cell  creates  a 
special  acoustic  environment  depending  on  its  geometrical  design.  Therefore,  a  second  objective 
of  this  primer  is  to  provide  a  simple  mathematical  analysis  of  the  acoustic  characteristics  of  a  test 
cell.  The  noise  source  and  the  acoustic  property  of  a  test  cell  may  under  certain  circumstances  be 
directly  coupled.  When  this  happens,  strong  resonance  may  occur.  This  form  of  interaction  tones 
is  also  discussed  in  this  primer. 

The  acoustic  spectrum  in  an  engine  test  cell  consists  of  tones  and  broadband  noise.  Whereas 
large  amplitude  tones  are  of  primary  concern  to  AEDC,  the  ability  to  predict  and  to  understand 
broadband  noise  as  well  as  medium  and  relatively  small  amplitude  tones  is  an  invaluable  asset  of 
a  test  engineer. 

The  source  of  acoustic  energy  in  an  engine  test  cell  is  the  jet  flow.  One  important  part  of  the 
noise  is  the  direct  noise  radiated  from  a  high-speed  turbulent  jet.  The  noise  components  are 
essentially  the  same  as  those  of  a  free  jet.  These  noise  components  are  discussed  in  Chapter  1. 
Another  important  noise  component  of  an  engine  test  cell  is  the  noise  generated  by  the  coupling 
of  the  jet  flow  to  the  acoustic  environment  provided  by  the  test  cell.  Usually,  the  coupling  will 
result  in  tones  occurring  at  the  resonance  frequencies.  This  type  of  resonant  tone  is  discussed  in 
Chapter  3.  Chapter  2  is  devoted  to  an  examination  of  the  acoustic  environment  existing  in  several 
types  of  engine  test  cells. 

In  Part  II  of  this  report,  a  method  to  suppress  super  resonance  is  proposed.  This  is  followed 
by  a  recommendation  on  the  establishment  of  a  computational  aeroacoustics  (CAA)  capability  at 
AEDC. 
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Part  I.  Aeroacoustics  of  Engine  Test  Cells 
Chapter  1.  The  Noise  of  High-Speed  Jets 

We  will  begin  by  performing  a  dimensional  analysis  of  jet  noise  experiments.  Sueh  an 
analysis  serves  two  useful  purposes.  First,  it  foeuses  attention  on  the  variables  or  physieal 
quantities  that  are  under  the  eontrol  of  an  experimentalist  or  a  test  cell  engineer.  In  addition,  such 
an  analysis  will  help  to  develop  scaling  formulas  so  that  measured  data  from  small  scale  models 
can  be  scaled  up  to  full  scale  test  cell  conditions.  Then,  we  will  turn  our  attention  to  the  three 
principal  components  of  supersonic  jet  noise.  Their  characteristics  will  be  examined  and 
prediction  methods  will  be  discussed. 


1.1  Dimensional  Analysis  of  Jet  Noise  Experiments. 

1.1.1  Derivation  of  Scaling  Formula 

We  will  consider  an  ideal  jet  noise  experiment  as  shown  in  Figure  1.1.  We  will  define  an 
ideal  experiment  as  one  that  satisfies  the  following  conditions/assumptions. 

1.  The  jet  velocity  profile  at  the  nozzle  exit  is  fairly  uniform  and  the  nozzle  internal  wall 
boundary  layer  is  reasonably  thin  such  that  it  has  no  strong  effect  on  the  radiated  jet  noise. 

2.  The  level  of  upstream  disturbances  is  normal  and  exerts  no  unusual  influence  on  the  noise  of 
the  jet. 

3.  Noise  absorption  by  ambient  air  (humidity  absorption)  can  be  ignored  or  the  measured  data 
has  been  adjusted.  In  other  words,  only  the  lossless  far  field  noise  is  considered. 

Conditions/assumptions  1  and  2  have  the  explicit  purpose  of  restricting  our  study  to  pure 
jet  mixing  noise.  Here,  to  avoid  misunderstanding,  we  wish  to  make  it  clear  that  an  ideal 
experiment  is  not  an  impossible  experiment.  Recently,  Viswanathan  &  Clark  (2004)  studied  the 
effects  of  internal  nozzle  boundary  layer  thickness  on  jet  noise.  They  concluded  that  unless  there 
was  flow  separation  inside  the  nozzle,  boundary  layer  thickness  had  minimal  impact  on  pure  jet 
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mixing  noise.  Condition/assumption  2  is  extremely  important.  It  has  been  known  sinee  the 
seventies  that  upstream  disturbanees  sueh  as  tones  eould  have  drastie  effeet  on  the  level  and 
speetral  shape  of  the  noise  of  a  jet.  In  a  more  reeent  study,  Viswanathan  (2003)  doeumented  the 
eontamination  of  jet  noise  by  various  eomponents  of  rig  noise  in  an  improperly  designed  test 
faeility.  But  onee  rig  noise  is  suffieiently  redueed,  the  quality  of  far  field  jet  noise  is  signifieantly 
improved,  opening  the  possibility  of  sealability.  Assumption  3  is  needed  beeause  humidity 
eorreetion  is  frequeney  and  henee  jet  size  speeifie.  Sueh  eorreetion  eannot  be  sealed. 

Poo  ,  Too  ,  Poo 


Figure  1.1.  Input  variables  of  an  ideal  jet  noise  experiment. 

We  will  refer  to  the  variables  of  an  experiment  that  affeet  the  noise  radiated  from  the  jet  as 
input  variables.  We  will  eall  variables  that  are  a  part  of  the  measured  or  proeessed  data  of  the 
experiment  “output  variables.”  The  far  field  noise  radiated  from  a  jet  is  infiueneed  by  both  the  jet 
exit  variables  as  well  as  the  variables  eharaeterizing  the  ambient  eonditions.  The  ambient 
variables  as  shown  in  Figure  1.1  are  poo,  and  Tk).  But  instead  of  Too,  we  may  use  aoo,  the 
ambient  sound  speed.  The  two  variables  are  equivalent  sinee  =  {yRToo)  .  Also  by  the  equation 
of  state  of  a  perfeet  gas,  we  have  =  ypj  ai,  so  that  only  two  of  the  three  ambient  variables  are 
independent.  We  will  take  pc^,  or  Ik  as  the  input  ambient  variables. 

The  jet  variables  are  Uj  (fully  expanded  jet  veloeity), (pressure  of  jet  at  the  nozzle  exit),  pj, 
Tj  (the  fully  expanded  jet  density  and  temperature),  v  (kinematie  viseosity)  and  Dj  (the  fully 
expanded  jet  diameter;  in  the  ease  of  subsonie  and  perfeetly  expanded  supersonie  jets,  it  is  the 
same  as  D,  the  nozzle  exit  diameter).  For  subsonic  and  perfectly  expanded  supersonic  jets,  the 
static  pressure  of  the  jet  at  the  nozzle  exit  is  almost  the  same  as  the  ambient  pressure.  That  is 
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Pe  =  Pco.  Thus  pe  is  not  an  independent  input  variable.  We  will  defer  the  eonsideration  of 
imperfectly  expanded  supersonic  jets,  for  which pe  is  not  equal  to  pcr>,  until  later. 

Again,  by  the  perfect  gas  law,  we  have  pj=  pJiRTj).  It  follows  that  pj  is  also  not  an 
independent  input  variable.  In  a  jet  noise  experiment,  it  is  more  convenient  to  control  the 
reservoir  temperature,  Tj-,  of  the  jet.  By  the  energy  equation,  Tj  is  related  to  Uj  and  Tj-  by. 


1  2 

C  T  =  C  T +—U 


Hence,  we  may  use  Tj-  instead  of  Tj  as  an  input  variable.  Therefore,  the  input  variables 
characterizing  the  jet  exit  conditions  are  Uj,  Tr,  uandZ);. 

Let  the  three  fundamental  dimensions  in  mechanics,  namely,  mass,  length  and  time  be 
denoted  by  M,  L  and  T,  respectively.  The  dimensions  of  all  the  input  variables  are  as  shown,  in 
parentheses,  below. 


( 

P\lT^)^ 


(£\ 


f 


T^RJ 


J 


Dp) 


where  R  is  the  gas  constant  of  the  equation  of  state.  It  is  to  be  noted  that  of  the  six  variables,  only 
Pec  has  dimension  M.  For  this  reason,  it  cannot  be  combined  with  other  input  variables  to  form  a 
dimensionless  group.  Three  independent  dimensionless  groups  can  be  formed  from  the 
remaining  five  variables.  A  convenient  choice  of  groupings  leads  to. 


t/,- 


V 


=  Re  (Reynolds  number  of  the  jet). 


Of  interest  in  a  jet  noise  experiment  is  the  far  field  noise  power  spectral  density  S'  at  a  point 
with  spherical  coordinates  (r,6(^zi)  (the  origin  of  the  coordinate  system  is  at  the  center  of  the 
nozzle  exit)  defined  by. 


/  =  j  S(r,0,<P,f)df . 
0 


(1.1) 
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where  is  the  time  average  of  the  square  of  the  pressure  fluetuations  at  (r,0,(l))  and  /is  the 
frequency.  A  simple  dimensional  analysis  of  (1)  indicates  that  the  dimensions  of  S  are 
At  the  same  time,  /  has  the  dimension  of  On  using  the  input  variables  to 
nondimensionalize  the  two  output  variables,  two  dimensionless  groups  can  be  formed.  They  are 


5 


and 


(Strouhal  number) 


Now  according  to  Buckingham  pi  theorem  (see  e.g.;  Sedov,  1959;  White,  2003),  the 
dimensionless  power  spectral  density  must  be  a  function  of  all  the  other  dimensionless  groups. 
Therefore,  we  may  write. 


SU^ 


F\ 


U: 


V^OO 


'TJ  f/,. 


.Re,6» 


kD.j 


(1.2) 


In  (1.2),  the  inverse  square  dependence  of  sound  on  the  distance  of  propagation  is  explicitly 
exhibited.  For  full  scale  jet  engines,  the  jet  Reynolds  number  is  in  the  millions.  At  such  a  large 
value,  the  function  F  would  be  insensitive  to  the  exact  value  of  the  Reynolds  number.  Another 
way  of  stating  this  is  that  F  would  have  essentially  reached  its  asymptotic  value  for  large  Re,  so 
that  it  is  practically  independent  of  Re.  This  reduces  (1.2)  to. 


F\ 


^U:  T  P  ^ 


’  F  ’  U, 


r 

v^y 


(1.3) 
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In  a  laboratory  experiment,  if  the  Reynolds  number  is  suffieiently  large,  say  half  a  million  or 
more,  (1.3)  would  also  apply.  When  this  is  true,  (1.3)  provides  a  basie  formula  for  sealing  small- 
seale  jet  noise  data  to  full  size  prototype. 

If  in  a  series  of  experiments,  the  ambient  pressure  is  praetieally  the  same,  then  we  may 
replaee  on  the  left  side  of  (1.3)  by  where  p^^^  is  the  reference  pressure  for  the  dB  scale 


plPi 


F 


,0 


r 


(1.4) 


However,  we  must  caution  that  (1.4)  is  just  an  approximation.  Eq.  (1.3)  is  the  correct  scaling 
formula  based  on  dimensional  analysis. 

1.1.2  Testing  Scaling  Formula 

Scaling  formulas  (1.3)  and  (1.4)  will  now  be  applied  to  a  variety  of  jet  noise  data.  The 
purpose  is  to  see  if,  indeed,  data  can  be  collapsed  according  to  these  formulas.  Also,  if  there  are 
problems  collapsing  the  data,  we  wish  to  examine  what  might  be  the  cause  of  the  problems. 

Figures  1.2(a)  and  1.2(b)  show  two  sets  of  data  from  Mach  2.0  jets  measured  by  Seiner  (see 
Tam,  Golebiowski  &  Seiner  1996).  One  jet  has  an  exit  diameter  of  3.60  in.  (9.14  cm),  the  other 
1.96  in.  (4.98  cm).  The  spectrum  data  are  plotted  according  to  formula  (1.4)  as  a  function  of 
Strouhal  number.  As  can  be  seen  from  these  figures,  there  is  a  good  collapse  of  data  in  all 
directions  of  radiation  in  spite  of  the  fact  that  the  smaller  jet  is  only  slightly  over  50%  of  the  size 
of  the  larger  jet.  At  Strouhal  number  above  6.0  the  spectra  of  the  smaller  jet  drops  off  abruptly. 
This  is,  most  likely,  an  instrumentation  problem.  This  part  of  the  spectra  should  be  disregarded. 

Figures  1.3(a)  and  1.3(b)  show  similar  comparisons  for  two  sets  of  Mach  0.5  data  at  a  jet 
temperature  ratio  of  2.15  measured  by  Bhat  (2001)  of  the  Boeing  company.  In  this  case,  the 
larger  jet  is  2.3  times  larger.  Over  the  angular  sector  from  6^=50°  to  6=  160°,  the  spectra 
collapse  fairly  well  except  at  very  low  and  very  high  Strouhal  numbers.  The  problem  at  low 
Strouhal  number,  where  the  spectra  of  the  large  jet  make  a  steep  drop  off,  is  due  primarily  to  the 
low  frequency  filter.  This  part  of  the  spectra  should  be  ignored.  At  high  Strouhal  number,  the 
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noise  spectra  of  the  smaller  jet  are  consistently  lower.  Figures  1.4(a)  and  1.4(b)  show  similar  data 
at  Mach  1.0.  Again  overall,  the  spectra  seem  to  scale  quite  close  to  each  other.  But  the  collapse 
again  is  not  good  at  very  low  and  very  high  Strouhal  numbers.  The  low  Strouhal  number 
problem  is  the  same  as  the  Mach  0.5  jet.  The  discrepancy  at  high  Strouhal  numbers  is  somewhat 
puzzling.  A  very  appealing  first  suggestion  is  that  it  is  a  Reynolds  number  effect.  The  Reynolds 
numbers  for  the  Mach  0.5  jets  are  1.87x10^  and  4.30x10^,  respectively.  The  smaller  jet  has  such 
a  low  Reynolds  number,  it  is  most  likely  that  the  initial  boundary  layer  is  transitional.  On  the 
other  hand,  the  larger  jet  having  a  Reynolds  number  close  to  half  a  million  should  be  nearly  fully 
turbulent.  One  would,  therefore,  expect  that  the  larger  jet  generates  more  high  frequency  noise. 
However,  for  the  Mach  1.0  jets,  the  Reynolds  numbers  are  4.36x10^  and  1.0x10^,  respectively. 
These  Reynolds  numbers  are  sufficiently  large  that  Reynolds  number  should  not  be  a  factor  in 
the  radiated  noise.  But  Figures  1.4(a)  and  1.4(b)  (as  well  as  Figure  1.2)  show  significant 
difference  in  the  noise  spectra  of  the  jets  at  high  Strouhal  number;  very  similar  to  those  of 
Figures  1.3(a)  and  1.3(b).  This  casts  doubt  on  the  proposition  that  Reynolds  number  is  really  the 
cause  of  the  discrepancy.  We  are,  at  the  present  time,  unable  to  find  a  good  explanation. 


Figure  1.2(a).  Scaling  of  jet  noise  spectra.  Data  from  Seiner.  Mj  =  2.0,  TrlT^^  =  1.8. 
- Dj  =  3.60  in.  (9.14  cm), . Dy  =  1.96  in.  (4.98  cm). 
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Figure  1.2(b). 


Figure  1.3(a).  Scaling  of  jet  noise  spectra.  Boeing  data  (Bhat  2001).  Mj  =  0.5, 
Tr/Too  =  2.15. - Dj  =  3.46  in.  (8.79  cm), . Dj=  1.5  m.  (3.81  cm). 
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Figure  1.4(a).  Scaling  of  jet  noise  spectra.  Boeing  data  (Bhat  2001).  Mj=  1.0, 
TJT^  =  2.15. - Dj  =  3.46  in.  (8.79  cm), . Z);  =  1.5  in.  (3.81  cm). 
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Recently,  Viswanathan  (2002)  measured  a  set  of  high  quality  jet  noise  data  over  a  wide  range 
of  Mach  numbers  and  temperature  ratios.  In  his  experiment,  nozzles  of  1.5,  2.45  and  3.46  inches 
or  3.81,  6.22,  and  8.79  cm  diameters  were  used.  Here  we  present  a  small  set  of  his  data  using 
(1.4)  as  the  scaling  formula.  Figures  5  and  6  are  noise  spectra  at  Tr/T^  =  1.8.  Mach  0.6  data  are 
shown  in  Figurel.5,  whereas  Mach  1.0  data  are  given  in  Figure  1.6.  Figures  1.7  and  1.8  show 
similar  spectra  at  a  much  higher  temperature  ratio  of  3.2.  It  is  evident  from  these  figures  that  the 
noise  spectra  from  the  two  larger  diameter  nozzles  do  collapse  well.  There  are  some  differences 
in  the  high  frequency  part  of  the  spectra  at  low  Mach  number  and  low  temperature  ratio  jet 
operating  conditions.  The  spectra  of  the  1.5  in.  (3.81  cm)  diameter  jet  are  consistently  lower  at 
high  Strouhal  number.  The  fact  that  the  spectra  of  the  two  larger  nozzles  collapse  into  a  nearly 
single  curve  suggests  that  the  collapsed  spectrum  should  be  very  close  to  the  pure  jet  noise 
spectrum  at  the  particular  jet  operating  conditions. 


statement  A:  Approved  for  public 
release;  distribution  is  unlimited. 


13 


AEDC-TR-06-3 


Figure  1.5(a).  Scaling  of  jet  noise  spectra.  Data  from  Viswanathan  (2002). 

Mj  =  0.6,  TrlT^=  1.8. - Dj  =  2.45  in.  (6.22  cm), . Dy  =  1.5  in. 

(3.81  cm). 


fD/U, 

Figure  1.5(b). 
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Figure  1.6(a).  Scaling  of  jet  noise  spectra.  Data  from  Viswanathan  (2002). 

Mj=  1.0,  Tr/Ta,=  \.S. - Dj  =  3.46  in.  (8.79  cm), - Dj  =  2.45  in.  (6.22 

cm), . Dy  =  1.5  in.  (3.81  cm). 
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Table  I 

Jet  Reynolds  number  at  =  3,2 


Mj\Dj 

1.5  inch  jet 

2.45  in. 

3.46  in. 

(3.81cm) 

(6.22  cm) 

(8.79  cm) 

0.6 

1.4x10^ 

2.38x10^ 

3.36x10^ 

1.0 

2.67x10^ 

4.50x10^ 

6.37x10^ 

Table  II 

Jet  Reynolds  number  at  TrITao  =  1,8 


Mj\Dj 

1.5  inch  jet 

2.45  in. 

3.46  in. 

(3.81  cm) 

(6.22  cm) 

(8.79  cm) 

0.6 

2.70x10^ 

4.59x10^ 

6.49x10^ 

1.0 

5.20x10^ 

8.80x10^ 

1.24x10^ 

The  reason  why  the  spectra  of  the  smallest  jet  of  Viswanathan’s  data  are  consistently  low  at 
high  Strouhal  number  is  not  known.  It  is  not  completely  a  Reynolds  number  problem.  Tables  I 
and  II  show  the  Reynolds  numbers  of  these  jets.  There  is  no  question  that,  as  shown  in  Table  I, 
the  Reynolds  number  is  low  at  temperature  ratio  3.2  for  the  smallest  jets.  However,  in  Table  II  at 
temperature  ratio  1.8  and  Mach  I.O,  the  Reynolds  number  of  the  smallest  jet  is  more  than  half  a 
million.  As  yet.  Figure  1.6  shows  that  the  noise  spectra  still  have  similar  discrepancies  at  high 
Strouhal  number  similar  to  the  lower  Reynolds  number  jet.  We  will  not  pursue  this  further  here, 
as  it  may  be  a  facility  problem  or  even  a  humidity  correction  problem.  Humidity  correction  starts 
to  become  important  at  20  KHz.  For  the  smallest  jet,  the  corresponding  Strouhal  number  lies  in 
the  range  of  1.3  to  2.7.  If  there  is  an  over  correction,  the  noise  spectra  would  be  lower  at  Strouhal 
number  higher  than  the  20  KHz  Strouhal  number.  We,  however,  have  no  evidence  that  this  is  the 
case. 
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1.1.3  Imperfectly  Expanded  Supersonic  Jets 


A 

Dj 

V 

Figure  1.9.  Schematic  diagram  of  an  underexpanded  supersonic  jet. 

For  imperfectly  expanded  supersonic  jets,  the  pressure  at  the  nozzle  exit  is  not  the  same  as 
the  ambient  pressure.  Thus  the  analysis  of  Section  11  is  not  directly  applicable  to  these  jets.  Let 
us  consider  an  underexpanded  jet  from  a  convergent-divergent  (C-D)  nozzle  as  shown  in  Figure 
1.9.  By  definition  the  static  pressure  of  the  jet  at  the  nozzle  exit  is  higher  than  the  ambient 
pressure.  The  gas  coming  out  of  the  nozzle  would,  therefore,  undergo  further  expansion.  This 
further  expansion  leads  to  an  effective  (or  fully  expanded)  jet  diameter,  Dj,  larger  than  the 
physical  diameter,  D,  of  the  nozzle  exit.  A  simple  relationship  between  Dj  and  D  can  be  found  by 
appealing  to  conservation  of  mass  flux. 

Since  the  gas  expands  downstream  of  the  nozzle  exit,  it  is,  therefore,  a  reasonably  good 
approximation  to  assume  that  the  flow  through  the  nozzle  as  well  as  in  the  initial  portion  of  the 
free  jet  is  isentropic.  We  will  use  subscripts  r,  e  and  j  to  denote  variables  in  the  reservoir,  at  the 
nozzle  exit  and  in  the  fully  developed  jet.  Conservation  of  enthalpy  leads  to, 

=  +  =  +  (1.5) 


Pa>  ,  Too  ,  Poo 


C-D  nozzle 


where  Cp  is  the  specific  heat  at  constant  pressure.  Eq.  (1.5)  may  be  rewritten  as. 


T. 


(pX" 

VPe) 


(1.6) 
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(1.7) 


The  last  equality  of  (1.6)  and  (1.7)  comes  from  isentropic  relation,  T  oc  f  In  (1.6)  Md  is  the 
design  Mach  number  of  the  CD  nozzle  {Md  =1.0  for  a  convergent  nozzle)  and  is  equal  to  the 
flow  Mach  number  at  the  nozzle  exit;  i.e.,  Md  =  Me.  Mj  in  (1.7)  is  the  fully  expanded  jet  Mach 
number.  The  combination  of  (1.6)  and  (1.7)  yields, 


K 

Pj 

and 


_j_ 

Te 


By  conservation  of  mass  flux,  we  have 

P,up=ppp.  (1.10) 

Thus,  by  using  (1.5),  (1.8)  to  (1.10),  we  find 


(1.9) 


l  +  ^Mj 


y-\ 


(1.8) 
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1  1 
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2(7-1) 


1  1 
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l  +  ^M] 


yd 


4(7*1) 


(1.11) 


Formula  (1.11)  was  first  stated  in  the  work  of  Tam  and  Tanna  (1982)  but  without  a  detailed 
derivation. 

We  want  to  point  out  that  for  an  overexpanded  jet,  the  jet  fluid  will  encounter  shocks 
immediately  downstream  of  the  nozzle  exit.  However,  for  moderately  overexpanded  jets,  the 
shocks  are  relatively  weak.  The  change  in  entropy  across  the  weak  shocks  of  these  jets  is  of 
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higher  order  in  the  jump  in  Maeh  number.  The  flow,  therefore,  does  not  deviate  badly  from 
isentropie  flow.  For  this  reason,  (1.11)  may  still  be  useful  for  estimating  the  fully  developed  jet 
diameter  of  moderately  overexpanded  jets. 

Shoek  eontaining  jets  usually  emit  sereeeh  tones  (Tam,  1995)  driven  by  a  feedbaek 
phenomenon.  For  strongly  sereeehing  jets,  the  flow  dynamies  are  substantially  ehanged  by  the 
feedbaek.  Here  we  eonsider  jets  that  have  only  mild  sereeeh  tones. 
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Figure  1.10.  Sealing  of  perfeetly  and  imperfeetly  expanded  supersonie  jet  noise 

speetra.  Mj=\.5.  -  Md=\.5,  Z)=1.68  in.(4.27  em),  . Mrf=1.0, 

D  =  1.56  in.(3.96  em)  TrITao  =  2.23  for  inlet  angle  90°,  TrlT,n  =  1.45  for  inlet  angle 
145°. 

Now  as  far  as  jet  mixing  noise  is  eoneemed,  Dj,  given  by  (1.11)  is  the  pertinent  length  seale. 
Thus  the  results  obtained  by  dimensional  analysis  before  would  apply  to  even  imperfeetly 
expanded  supersonie  jets  if  Dj  instead  of  D  is  used  in  the  sealing  formula  (1.3).  To  demonstrate 
the  applieability  of  sealing  formula  (1.3)  with  this  modifieation,  we  will  show  that  the  turbulent 


statement  A:  Approved  for  public 
release;  distribution  is  unlimited. 


21 


AEDC-TR-06-3 


mixing  noise  of  an  underexpanded  jet  ean  be  sealed  to  that  of  a  perfeetly  expanded  jet  by  means 
of  (1.3).  Figure  1.10  shows  two  sets  of  jet  noise  speetra  measured  by  Seiner  (see  Tam  et  al. 
1996)  at  fully  expanded  Maeh  number  1.5.  One  set  of  speetra  is  from  a  jet  issued  from  a 
eonvergent  nozzle  (M^=  1.0).  The  other  set  is  from  a  nearly  perfeetly  expanded  jet  {Md=  1.5). 
As  ean  be  seen  by  using  the  fully  expanded  jet  diameter  Dj  provided  by  (1.11),  there  is  a  good 
eollapse  of  the  data.  The  underexpanded  jet  has  additional  broadband  shook  assooiated  noise  at 
high  Strouhal  number  at  0  =  90°  but  not  at  0  =  145°.  Figure  1.11  shows  similar  data  sealing  for 
Maeh  2  jets.  In  this  ease,  Dj  is  muoh  larger  than  D  for  underexpanded  jets.  The  speetra  eollapse 
well  using  Dj  of  (1.11)  as  the  length  soale.  The  data  level  will  differ  by  several  dB  if  D  is  used 
instead. 
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Figure  1.11.  Sealing  of  perfeetly  and  imperfeotly  expanded  supersonio  jet  noise 

speetra.  Mj  =  2.0  Tr/Too=  1.8.  - Md  =  2.0,  D=  1.96  in.(4.99  om), . 

Md=  1.0,D=  1.56in.(3.96  om) 
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1.2  Power  Law 

It  has  been  known  sinee  the  early  days  of  jet  noise  research  that  the  noise  intensity  varies  as  a 
high  power  of  the  jet  velocity.  This  power  law  has  the  support  of  both  theory  and  experiment. 
The  most  celebrated  power  law  is  the  Lighthill  if  law.  Numerous  experiments  have  been 
performed  for  the  purposes  of  verifying  the  if  law.  Some  measurements  indicate  that  the 
exponent  is  close  to  but  not  precisely  equal  to  8.  A  more  general  power  law  is  sometimes  written 
in  the  form, 


/oc[/;.  (1.12) 

In  a  recent  experiment,  Viswanathan  (2002)  confirmed  that  jet  noise  intensity  could,  indeed,  be 
correlated  by  a  power  law  of  the  form  (1.12).  He  found  that  n  depends  on  jet  temperature.  It 
varied  from  7.98  to  8.74. 

A  power  law  in  the  form  of  (1.12)  is  dimensionally  unsatisfactory.  In  this  section,  our  goal  is 
to  express  the  power  law  in  a  dimensionless  form  and  then  examine  the  dependence  of  its 
parameters  on  the  dimensionless  groups  of  the  input  variables. 

Let  us  define  the  noise  intensity  (or  OASPL,  overall  sound  pressure  level),  /,  for  noise 
radiation  in  the  far  field  direction  6  as, 

CO 

I  =  \s{rAf)if.  (1.13) 

0 

Note;  noise  intensity  defined  by  (1.13)  differs  from  the  traditional  definition  by  a  factor  of  p^aa- 
In  (1.13)  5*  is  the  spectral  density.  Now  (1.13)  may  be  rewritten  in  a  dimensionless  form  as. 
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On  replacing  the  integrand  by  (1.3),  (1.14)  becomes. 
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Experimental  observations  suggest  that  the  funetion  K  may  be  expressed  as  a  high  power  of 
(Uj/aoD).  That  is, 


pI 


A 


kD,j 


(1.16) 


where  A  and  n  are  two  free  parameters  of  the  power  law.  They  depend  on  the  temperature  ratio 
and  the  direetion  of  radiation;  i.e.,  A  =  A{Tr/Tao,0),  n  =  n{TrlT^,0). 

The  validity  of  power  law  in  the  form  of  (1.16)  ean  be  eonfirmed  experimentally.  Figure  1.12 
shows  a  log-log  plot  of  // pi,  at  rID  =  100.0  as  a  funetion  of  Ujla«,  for  jet  noise  radiated  in  the 
90°  direetion  based  on  the  data  of  Viswanathan  (2002)  and  Seiner  (see  Tam  et  al.  1996)).  Eaeh 
of  the  data  sets  is  for  a  fixed  temperature  ratio.  The  arrow  marks  the  value  -80.0  in  the  vertieal 
axis  of  Figure  1.12  for  the  partieular  eurve.  It  is  elear  that,  to  a  good  degree  of  approximation,  a 
linear  relationship  exists  in  the  log-log  plot.  This  validates  a  power  law  relationship  as  given  in 
(1.16).  Figure  1.13  shows  a  similar  plot  for  noise  radiated  in  the  150°  inlet  angle  direetion.  Again 
a  linear  relationship  between  I ! pi,  and  Uj/aa:,  exists  at  eaeh  temperature  ratio.  In  addition  to  the 
data  shown  in  Figures  1. 1 2  and  1. 1 3,  we  would  like  to  report  that  we  have  examined  data  at  other 
angles  (from  50°  to  150°).  All  of  them  ean  be  eorrelated  by  a  power  law.  We  believe  this 
provides  strong  empirieal  support  for  the  general  validity  of  (1.16). 

The  exponent  n  of  power  law  (1.16),  whieh  is  the  slope  of  the  straight  lines  in  Figures  1.12 
and  1.13,  varies  slightly  with  jet  temperature  ratio  and  the  direetion  of  radiation.  The  same  is  true 
with  the  eoeffieient  A.  Figure  1.14  gives  the  dependenee  of  n  on  temperature  ratio  at  a 

number  of  direetions  of  radiation.  Generally  speaking,  n  is  large  for  larger  inlet  angle.  Its  value 
deereases  with  inerease  in  temperature  ratio.  Numerieally,  n  varies  from  5.3  to  9.9  and  u  =  8  may 
be  eonsidered  an  overall  averaged  value. 

The  dependenee  of  eoeffieient  A  of  power  law  (1.16)  on  jet  temperature  ratio  is  given  in 
Figure  1.15.  Just  as  for  exponent  n,  the  value  of  A  is  larger  for  larger  inlet  angle.  Also  its  value 
deereases  nearly  monotonieally  with  inerease  in  temperature  ratio. 
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Figure  1.12.  Dependence  of  noise  intensity  on  jet  Mach  number  based  on  ambient 
speed  of  sound.  Inlet  angle  90°.  <—  vertical  axis  equals  to  -80.0  for  the  curve.  #, 

V,  O,  O,  A;  A/Foo  =  1.0,  1.8,  2.2,  2.7,  3.2. - Power  law  A  {Ujla^fl{rlDf. 

Data  from  Viswanathan  (2002)  and  Seiner  (see  Tam  et  al.,  1996). 


U/a 


Figure  1.13.  Dependence  of  noise  intensity  on  jet  Mach  number  based  on  ambient 
speed  of  sound.  Inlet  angle  150°.  <—  vertical  axis  equals  to  -80.0  for  the  curve.  #, 

V,  O,  O,  A;  A/Too  =  1.0,  1.8,  2.2,  2.7,  3.2. - Power  law  A  {Ujla^fl{rlDf. 

Data  from  Viswanathan  (2002)  and  Seiner  (see  Tam  et  al.  1996). 
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It  is  interesting  to  point  out  that,  as  shown  in  Figures  1.14  and  1.15,  n  and  A  deerease 
monotonieally  as  temperature  ratio,  TrlT^^,  inereases.  However,  both  n  andH  eannot  beeome  zero. 
Thus,  it  is  reasonable  to  expeet  that,  as  temperature  ratio  beeomes  suffieiently  large,  n  and  A 
would  eaeh  independently  reaeh  an  asymptotie  value.  Figures  1.14  and  1.15  suggest  that  at  the 
last  data  point  {TylT^  =  3.2)  both  eurves  for  n  and  A  are  elose  to  their  asymptotes.  We  will  eall 
this  the  “hot-jet”  limit.  The  physieal  signifieanee  of  the  hot-jet  limit  is  that  for  jets  at  high 
temperature  ratio,  the  noise  intensity  depends  only  on  jet  veloeity  ratio.  Temperature  plays  no 
role  in  jet  noise  intensity  exeept  through  its  influenee  on  the  nozzle  exit  veloeity.  This  is  an 
unexpeeted  result.  Of  eourse,  it  would  require  further  investigation  before  one  ean  determine 
what  the  hot-jet  limit  means  to  the  distribution  of  jet  noise  sourees  in  the  jet  plume  and  on  the  jet 
mixing  proeesses. 

1,3  The  Principal  Components  of  Supersonic  Jet  Noise 


St 

Figure  1.16.  Typieal  far  field  narrow  band  supersonie  jet  noise  speetrum  (Seiner 
1984).  Mierophone  at  30  degrees  to  nozzle  inlet  direetion. 
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Except  for  jets  operating  at  perfeetly  expanded  eonditions,  the  noise  of  a  supersonie  jet 
eomprises  three  basie  eomponents.  They  are  the  turbulent  mixing  noise,  the  broadband  shook 
assooiated  noise  and  the  soreeoh  tones.  Figure  1.16  shows  a  typioal  narrow  band  noise  speotrum 
of  an  imperfeotly  expanded  supersonie  jet  measured  by  Seiner  (1984).  Very  prominent  in  this 
figure  is  the  disorete  frequenoy  oomponent  or  the  soreeoh  tone  at  the  oenter  of  the  figure.  The 
appearanoe  of  a  soreeoh  tone  is  usually  aooompanied  by  its  harmonies.  Sometimes,  even  the 
fourth  or  fifth  harmonio  oan  be  deteoted.  The  low  frequenoy  broadband  peak  of  Figure  1.16,  to 
the  left  of  the  soreeoh  tone,  is  the  turbulenoe  mixing  noise  oomponent.  The  high  frequenoy  peak, 
to  the  right  of  the  fundamental  soreeoh  frequenoy,  is  the  broadband  shook  assooiated  noise.  The 
relative  intensity  of  the  three  noise  eomponents  is  a  strong  funotion  of  the  direotion  of 
observation.  In  the  downstream  direotion  of  the  jet,  turbulent  mixing  noise  is  the  most  dominant 
noise  oomponent.  In  the  upstream  direotion,  the  broadband  shook  assooiated  noise  is  more 
intense.  For  oiroular  jets,  the  soreeoh  tones  radiate  primarily  in  the  upstream  direotion. 

1.3.1  Characteristics  of  Turbulent  Mixing  Noise 

High  quality  turbulent  mixing  noise  data  in  narrowband  from  nearly  perfeetly  expanded 
supersonie  jets  oan  be  found  in  the  literature.  Seiner  et  al.  (1992)  performed  very  detailed  narrow 
band  noise  measurements  of  a  Maoh  2  jet.  The  total  temperature  of  the  jet  varied  from  313K  to 
as  high  as  1534K.  Figure  1.17  shows  typioal  noise  direotivities  at  a  few  of  the  important  Strouhal 
numbers  (Strouhal  number  =fDIUj,  where  f,  D  and  Uj  are  the  frequenoy,  nozzle  exit  diameter 
and  the  fully  expanded  jet  velooity,  respeotively).  It  is  olear  from  this  figure  and  all  the  other  data 
that  the  dominant  part  of  the  turbulent  mixing  noise  is  radiated  into  an  angular  seotor  of  about  45 
degrees  measured  from  the  jet  flow  direotion  or  greater  than  135  degrees  inlet  angle.  Within  this 
dominant  noise  seotor,  the  noise  speotrum  is  dominated  by  a  single  broadband  peak.  The  peak 
ooours  at  a  Strouhal  number  of  around  0.1  to  0.25  depending  on  the  jet  temperature  and  Maoh 
number.  Outside  the  dominant  noise  radiation  seotor,  the  noise  direotivity  and  speotral 
oharaoteristios  are  very  different.  In  the  upstream  direotion,  the  noise  intensity  is  low  and  is 
nearly  uniform  in  direotion  as  shown  in  Figure  1.17.  Furthermore,  the  noise  speotrum  is  very  fiat 
over  a  large  range  of  Strouhal  number.  This  is  the  baokground  noise,  whioh  is  believed  to  be 
generated  by  the  fine  soale  turbulenoe  of  the  jet  flow.  The  dominant  part  of  the  noise,  on  the 
other  hand,  is  believed  to  be  generated  direotly  by  the  large  turbulenoe  struotures  of  the  jet  flow. 
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X,  deg 

Figure  1.17.  Measured  noise  direetivities  (Seiner  et  al.,  1992)  at  seleeted  Strouhal 
numbers  (St)  of  a  Maeh  2  jet  at  total  temperature  of  500K.  O  St  =  0.067;  □ 

St  =  0.12;  A  St  =  0.20;  O  St  =  0.40. 

The  intensity  as  well  as  the  direetional  and  speetral  eharaeteristies  of  the  turbulent  mixing 
noise  of  supersonie  jets  depends  on  the  jet  Maeh  number  and  jet  to  ambient  temperature  ratio.  As 
jet  Maeh  number  inereases,  the  differenee  in  peak  noise  levels  between  the  dominant  part  and  the 
baekground  noise  inereases.  However,  at  a  fixed  Maeh  number  and  Strouhal  number,  the 
maximum  level  of  noise  radiation  does  not  signilioantly  inerease  with  large  inerease  in  jet 
temperature.  The  main  effeet  of  higher  jet  temperature  is  the  broadening  of  the  dominant  noise 
seetor. 

1.3.2  Characteristics  of  Broadband  Shock  Associated  Noise 

When  a  supersonie  jet  is  operating  at  an  imperfeetly  expanded  eondition,  regardless  of 
whether  the  jet  is  underexpanded  or  overexpanded,  a  shoek  eell  strueture  will  form  in  the  jet 
plume.  In  the  presenee  of  the  shoek  eells,  the  jet  emits  additional  broadband  noise  ealled  the 


statement  A:  Approved  for  public 
release;  distribution  is  unlimited. 


29 


AEDC-TR-06-3 


broadband  shock  associated  noise.  Broadband  shook  assooiated  noise  was  first  identified  by 
Harper-Bourne  &  Fisher  (1974).  Sinoe  then,  it  has  been  the  subjeot  of  a  number  of  experimental 
investigations;  e.g.,  Tanna  (1977),  Seiner  &  Norum  (1979,  1980),  Seiner  &  Yu  (1982),  Norum  & 
Seiner  (1982a,b)  and  Yamamoto  et  al.  (1984).  These  experiments  have  provided  muoh  detailed 
information  on  the  oharaoteristios  of  broadband  shook  assooiated  noise. 


0  5  10  15  20  25  30 

f  (KHz) 

Figure  1.18.  Narrow  band  noise  speotra  of  a  Maoh  1.67  jet  from  a  C-D  nozzle 
with  design  Maoh  number  1.5.  x°=  nozzle  inlet  angle  (Norum  &  Seiner  1982a). 

Figure  1.18  shows  typioal  noise  speotra  of  an  imperfeotly  expanded  jet  measured  by  a 
miorophone  array  at  equal  distanoe  to  the  jet  nozzle  exit  at  15°  intervals.  The  maximum  intensity 
of  eaoh  speotrum  is  indioated  in  the  figure.  The  jet  emitted  a  soreeoh  tone  at  a  frequenoy  of 
3.6KHz.  The  soreeoh  tone  oan  be  seen  in  almost  all  the  speotra.  The  turbulent  mixing  noise 
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contributes  to  the  broadband  peak  to  the  left  of  the  sereeeh  frequeney.  The  very  distinet  peak  to 
the  right  of  the  sereeeh  tone  is  the  broadband  shook  assooiated  noise.  Clearly,  this  is  the 
dominant  noise  oomponent  in  the  upstream  direotion.  Several  distinet  oharaoteristios  of 
broadband  shook  assooiated  noise  oan  be  identified  in  this  figure.  First,  the  frequeney 
oorresponding  to  the  peak  of  the  speotrum  ohanges  with  the  direotion  of  radiation.  It  inoreases 
with  inoreasing  y**,  the  inlet  angle.  Seoond,  eaoh  speotrum  is  made  up  of  several  peaks.  A  seoond 
and  a  third  peak  oan  readily  be  identified  in  the  ~  90  degree  speotrum.  Third,  the  half-width  of 
the  speotral  peak  inoreases  with  y*’.  Thus  at  large  the  speotral  peak  is  very  broad.  One  very 
unexpeoted  oharaoteristio  is  that  at  a  fixed  jet  Maoh  number,  the  maximum  noise  level  at  the 
speotral  peak  is  unaffeoted  by  jet  temperature.  The  frequeney  at  whioh  the  peak  level  ooours  is, 
however,  shifted  to  a  higher  value  as  jet  temperature  inoreases. 


1.3.3  Characteristics  of  Screech  Tones 


Mj 

Figure  1.19.  Sereeeh  modes  of  axisymmetrio  supersonio  jet  from  a  C-D  nozzle 
with  design  Maoh  number  1.41  (Seiner  eta/.  1987). 


Sereeeh  tones  from  supersonio  jets  were  first  observed  by  Powell  (1953a,b).  Sinoe  then  the 
phenomenon  has  been  studied  experimentally  by  a  host  of  investigators;  e.g.,  Davies  &  Oldfield 
(1962),  Sherman  et  al.  (1976),  Westley  &  Wooley  (1969,  1970,  1975),  Rosfjord  &  Toms  (1975), 
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Norum  (1983),  Seiner  et  al.  (1987),  Massey  et  al.  (1994)  and  more  reeently  Panda  (1998),  Panda, 
Seasholtz  (1999),  Panda,  Raman  &  Zaman  (2003).  It  is  found  that  the  fundamental  sereeeh  tone 
radiates  primarily  in  the  upstream  direetion  whereas  the  prineipal  direetion  of  radiation  of  the 
seeond  harmonie  is  at  90  degrees  to  the  jet  flow  direetion  (Norum,  1983).  When  a  jet  emits 
strong  sereeeh  tones  the  jet  flow  undergoes  strong  oseillations.  Observations  reveal  that  two 
types  of  jet  motion  eould  oeeur.  They  are  the  toroidal  mode  and  the  helieal/flapping  mode 
oseillations.  (Note;  the  superposition  of  an  equal  number  of  left  and  right  helieal  mode 
oseillations  ereates  a  flapping  mode.)  With  the  toroidal  mode  oseillations,  both  the  jet  flow  and 
the  aeoustie  field  exhibit  axisymmetry.  With  the  flapping  mode,  the  jet  oseillates  up  and  down 
around  a  flapping  plane.  The  flapping  plane  may  not  neeessarily  be  stationary.  A  slow  preeession 
of  the  flapping  plane  has  been  observed. 

Figure  1.19  shows  the  variation  of  the  aeoustie  wavelength  of  the  sereeeh  tones  of  a 
supersonie  jet  issued  from  a  Maeh  1.41  eonvergent-divergent  nozzle  (Seiner  et  al.  1987)  as  a 
funetion  of  jet  Maeh  number.  Four  modes  of  sereeeh,  two  toroidal  modes  and  two 
helieal/flapping  modes,  were  observed.  At  eertain  jet  Maeh  numbers,  the  jet  emits  more  than  one 
sereeeh  tone  (not  harmonie  of  eaeh  other).  The  dominant  sereeeh  tone  mode  ehanges  with  jet 
Maeh  number.  This  mode  switehing  phenomenon  is  elearly  shown  in  Figure  1.19.  At  low  Maeh 
number,  the  sereeeh  tones  are  assoeiated  with  the  toroidal  mode.  As  Maeh  number  inereases, 
there  is  a  switeh  over  to  the  flapping/helieal  mode. 

The  intensities  of  sereeeh  tones  are  affeeted  by  the  jet  Maeh  number,  the  jet  temperature,  the 
nozzle  lip  thiekness  and  the  presenee  of  sound  refleeting  surfaees  in  the  immediate  environment 
of  the  jet.  It  is  well-known  that  the  intensity  of  sereeeh  tones  deereases  with  inerease  in  jet 
temperature.  This  is  espeeially  true  for  low  supersonie  Maeh  number  jets  (see  Massey  et  al. 
1994).  Perhaps  beeause  of  this,  engine  designers,  generally,  have  not  regarded  sereeeh  tones  as  a 
potential  eause  of  sonie  fatigue  to  jet  engines.  A  thiek  nozzle  lip  ean  aet  as  a  good  sound 
refleeting  surfaee.  Norum  (1983)  found  that  the  sereeeh  intensity  ean  inerease  substantially 
simply  by  using  a  thieker  lip  nozzle.  The  tone  direetivity  is,  however,  not  altered. 

1,4  Turbulent  Mixing  Noise 

There  is  now  a  growing  aeoeptanee  that  turbulent  mixing  noise  of  high-speed  jets  eonsists  of 
two  eomponents.  They  are  the  noise  from  the  fine  seale  turbulenee  and  the  large  turbulenee 
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structures  of  the  jet  flow.  Laufer  et  al.  (1976)  used  a  spherical  reflector  to  measure  the  noise 
source  along  the  length  of  supersonie  jets  at  Maeh  1.47,  1.97  and  2.47.  They  found  that  the  noise 
source  strength  distributions  of  these  jets  radiating  to  the  maximum  noise  direetion  (around  30° 
to  45°  exhaust  angles)  and  to  the  sideline  at  90  degrees  were  distinetly  different.  Real  time 
pressure  signal  measured  by  an  omni  directional  microphone  indieated  that  the  sound  radiated  to 
those  two  directions  were  also  distinctly  different.  The  pressure  signal  at  90  degrees  resembled 
that  of  subsonie  jets  being  random  but  smooth.  On  the  other  hand,  for  the  two  high  Mach  number 
jets,  the  pressure  signal  at  37.5  and  45  degrees  had  shoek-like  spikes  with  steep  rise  time.  These 
distinctive  differences,  together  with  the  difference  in  source  locations,  led  Laufer  et  al.  to 
eonclude  that  there  were  two  intrinsieally  different  noise  sources.  The  noise  radiated  by  the  two 
sources  had  distinetively  different  direetional,  speetral  and  real  time  charaeteristies. 

Independently,  Tam  (1995),  Tam  and  Chen  (1993),  based  on  the  observation  of  a  strong 
direetivity  of  the  radiated  sound  in  the  low  exhaust  angle  direetions  (within  a  eone  of  50  degrees 
around  the  jet  flow  direction)  and  an  almost  flat  directivity  for  sound  radiating  upstream  and  to 
the  sideline  (see  Figure  1.17),  together  with  theoretical  reasoning  of  possible  noise  generation 
mechanisms,  proposed  that  the  dominant  part  of  the  noise  radiated  downstream  was  from  the 
large  turbulenee  structures  while  the  sideline  noise  was  from  the  fine  scale  turbulenee  of  the  jet 
flow.  During  the  last  few  years,  experimental  supports  are  found  validating  the  two  mixing  noise 
concepts. 

1.4.1  Noise  Generation  Mechanisms 

The  fine  seale  turbulenee  in  a  jet  flow  is  quite  isotropie.  The  random  turbulent  motion  creates 
an  effective  pressure  analogous  to  the  pressure  created  by  the  random  molecular  motion  of  a  gas. 
The  random  fluctuating  pressure  is  the  source  of  fine  scale  turbulence  noise.  Just  as  in  a  gas,  the 
turbulence  pressure  is  quite  isotropic.  As  a  result,  the  radiated  sound  has  a  nearly  uniform 
directivity.  The  directivity  is,  however,  modified  by  the  downstream  convection  of  the  fine  scale 
turbulenee  by  the  jet  mean  flow.  This  eonvection  effeet  eauses  more  sound  to  radiate  in  the 
direction  of  motion.  In  addition,  the  sound  generated  by  a  lump  of  fine  scale  turbulenee  has  to 
traverse  the  jet  mixing  layer  before  radiating  to  the  far  field.  Beeause  the  jet  mean  flow  is  non- 
uniform,  the  acoustie  rays  are  bent  due  to  refraction.  This  reduces  the  intensity  of  sound  radiated 
downstream. 
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Figure  1.20.  Schematic  diagram  showing  Mach  wave  radiation  generated  by  a 
wavy  surface  traveling  at  supersonic  phase  velocity. 

Let  us  now  examine  how  sound  is  generated  by  the  large  turbulence  structures  of  supersonic 
jets.  The  large  turbulence  structures  may  be  regarded  as  a  superposition  of  a  spectrum  of 
instability  waves  of  the  jet  flow.  The  two  are  statistically  equivalent.  Consider  an  instability 
wave  of  a  certain  azimuthal  mode  at  a  given  frequency  propagating  down  a  jet  column.  An 
approximate  picture  of  the  physical  problem  is  to  regard  the  instability  wave  as  a  wavy  wall.  The 
wavy  wall  has  the  same  wavelength  and  wave  speed  as  the  instability  wave.  In  the  case  of  the 
wavy  wall,  it  is  well-known  that  if  the  propagation  speed  is  supersonic  relative  to  the  ambient 
sound  speed,  intense  noise  radiation  in  the  form  of  Mach  waves  would  result  as  illustrated  in 
Figure  1.20.  The  wavy  wall  analogy  suggests  that  the  direction  of  the  most  intense  noise 
radiation  from  a  supersonic  jet  can  be  estimated  by  using  the  Mach  angle  relation  based  on  the 
speed  of  the  most  amplified  instability  wave  of  the  jet.  Further,  the  Strouhal  number  at  the  peak 
of  the  radiated  noise  spectrum  must  be  the  same  as  the  Strouhal  number  of  the  most  amplified 
instability  wave.  Tam  et  al.  (1992)  investigated  these  relationships  and  found  good  experimental 
confirmation. 

Tam  &  Burton  (1984)  pointed  out  that  the  wavy  wall  analogy  must  be  modified  to  account 
for  the  growth  and  decay  of  the  instability  wave  as  it  propagated  downstream.  Near  the  nozzle 
exit,  the  mixing  layer  of  the  jet  is  thin.  The  large  velocity  gradient  across  the  shear  layer  causes 
the  instability  wave  to  have  very  large  growth  rate.  As  the  instability  wave  propagates 
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downstream  the  shear  layer  thiekness  inereases  leading  to  a  reduetion  in  mean  veloeity  gradient 
and  henee  the  growth  rate.  Eventually  at  some  point  downstream  the  wave  beeomes  damped.  Its 
amplitude  deereases  as  it  eontinues  to  propagate  until  it  beeomes  vanishingly  small.  The  growth 
and  deeay  of  the  wave  amplitude  is  important  to  the  noise  radiation  proeess.  For  a  fixed 
frequeney  wave  of  eonstant  amplitude,  the  wave  speetrum  is  diserete.  With  a  single  wavenumber 
there  is  only  a  single  wave  speed  so  that  the  Maeh  waves  are  radiated  in  a  single  direetion.  The 
growth  and  deeay  of  the  instability  wave  amplitude  lead  to  a  broadband  wavenumber  speetrum. 
This  results  in  Maeh  wave  radiation  over  large  angular  direetions.  Furthermore,  a  single 
frequeney  subsonie  wave  of  eonstant  amplitude  would  not  radiate  sound  aeeording  to  the  Maeh 
wave  radiation  meehanism.  However,  with  growth  and  deeay  of  the  wave  amplitude,  a  part  of  the 
broadband  wavenumber  speetrum  eould  have  supersonie  phase  veloeity.  These  supersonie  phase 
disturbanees  will  lead  to  noise  radiation. 


Figure  1.21.  Sehematie  diagram  of  the  direetivity  of  Maeh  wave  radiation  and  that 
of  the  fine  seale  turbulenee. 


The  Maeh  wave  radiation  meehanism  diseussed  above  relies  on  the  existenee  of  supersonie 
phase  (relative  to  ambient  sound  speed)  eomponents.  For  highly  supersonie  jets,  espeeially  at 
high  temperature,  this  is  an  extremely  effieient  noise  generation  proeess.  But  if  the  jet  speed  is 
subsonie,  the  effieieney  is  greatly  redueed.  Thus  for  low  subsonie  jets,  the  large  turbulenee 
struetures  of  the  jet  flow  may  not  be  the  most  important  noise  generators.  The  fine  seale 
turbulenee  probably  is  the  more  dominant  souree. 
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Figure  1.21  is  a  schematic  diagram  showing  the  radiation  of  the  two  mixing  noise 
components  from  a  high  speed  jet.  Within  the  Mach  cone  in  the  downstream  direction,  the  large 
turbulence  structures’  noise  is  dominant.  In  the  upstream  and  sideline  directions,  there  is  little 
Mach  wave  radiation.  The  dominant  noise  is  from  the  fine  scale  turbulence. 

1.4.2  Experimental  Evidence  in  Support  of  Two  Turbulent  Mixing  Noise  Components 


Figure  1.22.  Similarity  spectra  for  the  two  components  of  turbulent  mixing  noise. 

- large  turbulence  structures/instability  waves  noise;  — - fine  scale 

turbulence  noise. 

Other  than  the  early  work  of  Laufer  et  al.  (1976),  an  important  set  of  experimental  evidence 
in  support  of  the  existence  of  two  turbulent  mixing  noise  components  in  high  speed  jets  was 
offered  by  Tam,  Golebiowski  and  Seiner  (1996).  They  analyzed  a  very  large  set  of  jet  noise  data 
from  the  Jet  Noise  Laboratory  of  the  NASA  Langley  Research  Center.  They  reasoned  that  jet 
flow  and  turbulence  were  known  to  exhibit  self-similarity.  It  was,  therefore,  possible  that  jet 
noise  spectra  also  possessed  similarity  characteristics.  They  identified  empirically  two  similarity 
spectra  that  seemingly  were  capable  of  fitting  all  noise  spectra  in  all  directions  of  radiation.  The 
two  spectra  F(/7/peak)  and  G(/7/peak)  are  shown  in  Figure  1.22.  The  is  the  frequency  at  the 
spectrum  peak.  The  F-spectrum  is  the  spectrum  of  the  large  turbulence  structures  noise.  It  fits  all 
jet  noise  spectra  radiated  within  the  Mach  cone  centered  at  the  jet  flow  direction  when  the  peak 
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of  the  spectrum  is  aligned  with  that  of  the  data.  This  is  true  regardless  of  the  direction  of 
radiation,  the  jet  Mach  number  and  temperature  ratio  of  the  jet.  This  is  demonstrated  in  Figure 
1.23.  In  this  figure,  four  spectra  from  two  different  Mach  numbers  at  several  different 
temperature  ratios  and  directions  of  radiation  are  shown  to  have  the  same  shape  as  the  F- 
spectrum. 


Figure  1.23.  Comparison  of  the  similarity  spectrum  of  large  turbulence 
structures/instability  waves  noise  and  measurements:  (a)  Mj  =  2.0,  r/Too  =  4.89, 

160.1°,  SPUax=  124.7  dB;  (b)My  =  2.0,  Tr/T^=l.l2,  ;r=  160.1°, 

SPLn,ax=  121.6  dB;  (c) My  =  1.96,  T,/Ta,=  l.7S,  z=  136.6°,  SPL^^ax  =  121.0  dB; 

(d)M,  =  1.49,  Tr/T^=\.\\,z=  138.6°,  SPLn^ax  =  106.5  dB. 

The  G-spectrum  is  the  spectrum  of  the  fine  scale  turbulence  noise.  It  fits  all  high-speed  jet 
noise  spectra  radiated  to  the  sideline  and  upstream  directions.  That  it  does  fit  the  measured 
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spectra  regardless  of  direction,  jet  Mach  number  and  temperature  ratio  can  be  seen  in  Figure 
1.24.  There  is  excellent  agreement  between  the  G-spectrum  and  the  measured  data.  In  the 
transition  direction  of  radiation  at  which  both  the  large  and  small  scale  turbulence  noise  are 
important,  it  is  necessary  to  use  a  combination  of  the  F  and  G  spectra  to  fit  the  measured  data. 
This  is  illustrated  in  Figure  1.25. 


Frequency,  Hz 


Figure  1.24.  Comparison  of  the  similarity  spectrum  of  large  turbulence 
structures/instability  waves  noise  and  measurements:  (a)  Mj=  1.49,  TrlT^  =  2.35, 
Z=92.9°,  SPUax  =  96.0  dB;  (b)M,  =  2.0,  T,/T^  =  4.89,  j=83.8°, 

SPLn,ax=  107.0  dB;  (c)  My  =  1.96,  Tr/T^  =  0.99,  ;r=  83.3°,  SPLmax  =  95.0  dB; 
(d)  My  =  1.96,  r/Too  =  0.98,  z=  120.2°,  SPLn^ax  =  100.0  dB. 
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Figure  1.25.  Comparison  of  the  sum  of  the  similarity  speetrum  of  both  large  and 
fine  seale  turbulenee  and  measurements  at  intermediate  direetion  of  radiation.  — 

- large  turbulenee  struetures  noise; - fine  seale  turbulenee  noise; - 

total:  (a) My  =  2.0,  TrlT^  =  Am,  160.1°,  SPhn^^x  =  124.7  dB;  (b) My  =  2.0, 
OToo=1.12,  ;^=  160.1°,  SPUax=  121.6  dB;  (o)My=1.96,  r,/roo=1.78, 

j=  136.6°,  SPLmax=  121.0  dB 

Aside  from  the  NASA  Langley  Researeh  Center  data,  Tam  (1998)  examined  the  noise 
speetra  of  supersonie  jets  issuing  from  simple  elliptie,  reetangular  and  plug  nozzles  and  found 
that  their  speetra  also  eonformed  to  the  similarity  speetra.  Tam  and  Zaman  (2000)  analyzed  the 
noise  speetra  of  eold  subsonie  jets  from  nonaxisymmetrie  and  tabbed  nozzles.  They  diseovered 
that  the  speetral  shapes  from  these  jets  had  the  same  shape  as  the  similarity  speetra.  Dahl  and 
Papamosehou  (2000)  reported  that  the  speetra  from  supersonie  eoaxial  jets  eould  also  be 
represented  by  the  two  similarity  speetra.  Most  reeently,  Viswanathan  (2000,  2004) 
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demonstrated  eonvineingly  that  these  two  speetra  also  were  good  fits  to  the  measured  speetra  of 
hot  subsonie  and  supersonie  jets  over  a  fairly  large  range  of  jet  operating  eonditions. 

Sinee  the  1960’s  there  have  been  many  attempts  to  measure  direetly  the  sourees  of  jet  mixing 
noise.  The  use  of  hot  wire  and  similar  instruments  led  to  flow  interferenee  and  has  doomed  many 
early  efforts.  Reeently,  Panda  et  al.  (2003)  sueeeeded  in  measuring  the  density  and  veloeity 
fluetuations  inside  a  jet  by  means  of  Rayleigh  seattering.  By  means  of  this  teehnique,  they  were 
able  to  measure  two-point  spaee-time  eorrelations  between  the  turbulent  fluetuations  at  a  point 
inside  the  jet  flow  and  the  aeoustie  pressure  measured  by  a  far  field  mierophone.  They  found  that 
when  the  mierophone  was  at  90  degrees  to  the  jet  flow  direetion,  there  was  very  little  eorrelation. 
However,  when  the  mierophone  was  moved  to  a  lower  exhaust  angle,  e.g.,  30  to  40  degrees, 
there  was  signifieant  eorrelation.  This  set  of  eorrelation  data  is  in  full  agreement  with  the 
existenee  of  two  noise  sourees,  the  large  and  fine  seale  turbulenee.  The  fine  seale  turbulenee  is 
random  and  spatially  uneorrelated.  Eaeh  souree  oeeupies  a  small  volume  produeing  a  small 
amount  of  sound.  Sinee  a  far  field  mierophone  reeeives  noise  from  all  sourees,  the  eorrelation  of 
its  pressure  signal  with  the  motion  of  a  blob  of  fine  seale  turbulenee  would  be  statistieally 
insignifieant.  This  is  the  ease  of  noise  radiating  to  the  sideline.  The  situation  is  totally  different 
for  a  mierophone  loeated  at,  say,  30°  from  the  jet  exhaust  direetion.  The  noise  in  this  direetion  is 
generated  by  the  large  turbulenee  struetures  of  the  jet.  In  other  words,  the  noise  souree  is 
spatially  eoherent  over  a  large  volume  of  the  jet  flow.  The  fluetuations  at  a  point  are  a  good 
representation  of  the  fluetuations  of  the  entire  eoherent  souree.  Thus,  it  is  no  mystery  that  there  is 
strong  eorrelation  of  the  single  souree  point  fluetuations  with  far  field  sound. 

1.4.3  Prediction  of  Fine  Scale  Turbulence  Noise  from  High-Speed  Jets 

An  aeeurate  predietion  theory  for  fine  seale  turbulenee  noise  from  high-speed  jets  has 
reeently  been  developed  by  Tam  and  Auriault  (1999).  The  Tam  and  Auriault  theory  is  semi- 
empirieal.  This  is  so  beeause  the  turbulenee  input  to  the  theory  is  not  from  first  prineiple.  There 
is  no  first  prineiple  turbulenee  theory  at  this  time.  In  their  theory,  the  souree  of  noise  is  the 
random  fluetuations  of  the  turbulent  kinetie  energy  of  fine  seale  turbulenee.  The  turbulenee 
intensity,  turbulenee  length  and  time  seales  are  given  by  the  k  -  s  turbulenee  model.  The  level  of 
empirieism  of  the  theory  is  at  about  the  same  level  as  the  k  -  ^  model. 
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The  original  Tam  and  Auriault  theory  was  developed  for  cold  to  moderate  temperature  jets.  It 
was  found  that  for  hot  jets,  the  standard  k-  ^  turbulence  model  is  inadequate.  Tam  and  Ganesan 
(2004)  proposed  a  modification  to  the  k  -  ^  model  to  include  the  effect  of  density  gradient  arising 
from  non-uniform  to  temperature  distribution.  This  improvement  of  the  k-  ^  model  allows  Tam, 
Pastouchenko  and  Viswanathan  (2005)  to  extend  the  Tam  and  Auriault  theory  to  very  high 
temperature  jets.  Let  x  denote  the  location  of  a  far  field  microphone.  With  respect  to  a  spherical 
polar  coordinate  system  centered  at  the  nozzle  exit  with  the  polar  axis  pointing  in  the  direction  of 
flow,  X  has  the  coordinates  {R,&,(/))  where  R  is  the  radial  distance,  and  0  and  (j)  are  the  polar  and 
azimuthal  angle.  The  radiated  noise  spectral  density  is  given  by. 
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where  Fjet  is  the  volume  of  the  jet  plume.  In  (1. 17),  r(v)  is  the  Gamma  function;  =  ^pk^,  ks  is 

the  intensity  of  the  kinetic  energy  of  fine  scale  turbulence  per  unit  volume;  fg  and  Ts  are  the 
corresponding  length  and  time  scales;  co  is  the  angular  frequency;  and  /7fl(x2,x,(2>)  is  the  adjoint 
Green’s  function  that  accounts  for  the  impedance  and  refraction  effects  of  the  jet  flow.  The  term 
c  is  an  empirical  constant.  The  spectral  density  S{\,a>)  is  per  unit  angular  frequency  co.  On 
converting  to  decibel  per  Strouhal  number  ifDjluj)  based  on  fully  expanded  jet  velocity,  uj,  and 
diameter,  Dj,  the  spectral  density  of  the  sound  field  at  {R,@,(p)  is 
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(1.18) 


where  />ref  is  the  reference  pressure  of  the  decibel  scale.  It  is  to  be  noted  that,  for  v=  A,  (I. 17) 
and  (1. 1 8)  reduced  to  the  spectral  density  formulas  of  Tam  and  Auriault. 
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The  noise  speetrum  formula  (1.17)  eontains  four  parameters,  namely,  q]lc^ ,  and  v.  For 
cold  to  moderate  temperature  jets,  explicit  dependence  of  these  parameters  on  the  k-s 
turbulence  model  is  given  in  Tam  and  Auriault.  For  hot  axisymmetric  jets,  the  density  gradient 
parameter  is  {l/p){dp/dr).  By  assuming  that  the  density  parameter  is  small,  its  first  order  effect 
may  then  be  represented  by  a  perturbation  term.  That  is,  it  may  be  represented  by  an  additional 
term  that  is  linear  in  {l/p){dp/dr).  On  balancing  the  dimensions  using  k  and  £  (k  and  £  are  the 
only  two  parameters  available  in  the  k-  ^  model  for  dimensional  adjustment),  Tam  et  al.  (2005) 
proposed  the  use  of  the  following  formulas. 
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(1.19(b)) 

(1.20(a)) 

(1.20(b)) 

(1.21(a)) 

(1.21(b)) 

(1.22(a)) 

(1.22(b)) 


In  (1.19)  to  (1.22)  the  (a)  formulas  are  to  be  used  when  du  ! dr  and  dp! dr  have  opposite  signs. 
The  (b)  formulas  are  to  be  used  when  they  have  the  same  sign.  The  reason  for  this  is  due  to  the 
effect  of  density  difference  on  the  Kelvin-Helmholtz  instability  in  the  shear  layer  of  the  jet.  Tam 
and  Ganesan  (2004)  have  pointed  out  that  enhanced  mixing  arising  from  density  effect  occurs 
only  when  the  lighter  fluid  is  moving.  The  values  of  k  and  £  in  these  formulas  are  found  by  the 
extended  k  -  £  turbulence  model  in  the  mean  flow  computation. 

For  nonaxisymmetric  jets,  (1.19)  to  (1.22)  may  be  generalized  by  replacing  \dpldr\  by 
|(VM)-(VyC>)  |/|  Vm  I  and  that  density  gradient  correction  is  incorporated  only  if  {Vu)-{Vp)  is 
negative.  Here  Vu  is  taken  as  the  reference  direction  as  far  as  density  gradient  effect  is 
concerned. 
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(14)  to  (17)  contained  four  unknown  constants.  They  are  c^,  B,  Cfp  and  c^p.  The 
recommended  values  of  these  constants  are, 


A  =  0.755,  C|  =  0.256,  =  0.233 

=  2.1599,  5  =  0.806, 

c,^  = -0.026,  c^^  = -0.2527 


Figure  1.26(a).  Cold  jet  noise  spectra  at  0  =  60°,  Tr/Too=l.0.  (a)  Mj  =  2.0, 

(b)  Mj  =  1 .49,  (c)  Mj  =  0.9,  (d)  Mj  =  0.7,  (e)  Mj  =  0.5,  (f)  Mj  =  0.3.  - 

experiments.  Seiner  (in  Tam,  Golebiowski  &  Seiner,  1996)  and  Norum  &  Brown 
(1993); . (1.17). 
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Figure  1.26(b).  Cold  jet  noise  spectra  at  0  =  90°,  TrlT^ 
experiments, . (1-17) 


Figure  1.26(c).  Cold  jet  noise  spectra  at  0=120°, 
experiments, . (1-17). 
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Figure  1.26  shows  comparisons  of  the  computed  fine  scale  turbulence  noise  spectra  by  Tam 
and  Auriault  (1999)  and  the  experimental  data  of  Seiner  (see  Tam,  Golebiowski  and  Seiner, 
1996)  and  Norum  and  Brown  (1993).  Figures  1.26(a),  1.26(b)  and  1.26(c)  show  the  noise  spectra 
at  60°,  90°  and  120°  for  cold  jets  with  Mach  number  ranging  from  0.3  to  2.0.  Notice  that  the 
frequency  range  spreads  over  3  decades.  The  sound  pressure  level  increases  by  60  dB  from  a 
Mach  0.3  jet  to  a  Mach  2.0  jet.  That  there  is  good  agreement  over  such  a  large  range  of 
frequency  and  sound  pressure  level  indicate  that  the  theory  must  contain  the  essential  physics  of 
the  generation  and  propagation  of  line  scale  turbulence  noise. 

Figures  1.27(a)  and  1.27(b)  show  comparisons  between  calculated  and  measured  noise 
spectra  at  inlet  angle  83°  and  93°  for  a  hot  Mach  2  jet.  The  data  are  from  the  NASA  Langley 
Research  Center  jet  noise  data  bank  of  Seiner.  The  spectra  at  three  temperatures  are  shown  in 
these  figures.  Temperature  ratio  TJTa  =  4.89  is  the  highest  known  temperature  ratio  at  which 
accurate  jet  noise  data  have  been  measured. 


Figure  1.27(a).  Comparisons  between  calculated  noise  spectra  and  experiment. 
Mj  =  2.0,  inlet  angle  =  83°.  (a)  T^/Ta  =  4.89,  (b)  VTa  =  4.08,  (c)  TJTa  =  3.28. 
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Figure  1.27(b).  Comparisons  between  calculated  noise  spectra  and  experiment. 
Mj  =  2.0,  inlet  angle  =  93°.  (a)  TJTa  =  4.89,  (b)  VTa  =  4.08,  (c)  VTa  =  3.28. 


1,5  Broadband  Shock  Associated  Noise 

Broadband  shock  associated  noise  and  screech  tones  are  generated  only  when  a  quasi- 
periodic  shock  cell  structure  is  present  in  the  jet  plume.  It  will  be  shown  that  the  quasi¬ 
periodicity  of  the  shock  cells  plays  a  crucial  role  in  defining  the  characteristics  of  both  the 
broadband  and  discrete  frequency  shock  noise.  The  importance  of  the  quasi-periodicity  of  the 
shock  cells  to  broadband  shock  associated  noise  was  first  recognized  by  Harper-Bourne  &  Fisher 
(1974)  in  their  pioneering  work  on  this  subject.  They  made  use  of  the  quasi-periodicity  to  explain 
some  of  the  prominent  characteristics  of  this  noise  component  discussed  earlier. 

1.5.1  Shock  Cell  Structure 

The  shock  cell  structure  of  an  imperfectly  expanded  supersonic  jet  is  formed  by 
oblique/normal  shocks  and  expansion  fans.  These  shocks  and  expansion  fans  are  generated  at  the 
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nozzle  lip  because  of  the  mismatch  of  the  static  pressures  inside  and  outside  the  jet.  For  an 
underexpanded  jet  an  expansion  fan  is  initiated  at  the  nozzle  lip.  The  fan  is  necessary,  for  it 
allows  the  static  pressure  to  decrease  gradually  to  that  of  the  ambient  condition.  For  an 
overexpanded  jet  an  oblique  shock  would  form  at  the  nozzle  lip  region.  In  passing  through  the 
shock  the  static  pressure  of  the  gas  increases  abruptly  to  match  that  outside  the  jet.  The 
expansion  fan  or  shock,  once  formed,  propagates  across  the  jet  flow  until  it  impinges  on  the 
mixing  layer  on  the  other  side.  Since  the  flow  outside  the  jet  is  stationary  or  subsonic  any  shock 
or  expansion  fan  is  not  allowed.  The  impinging  shock  or  expansion  fan  is,  therefore,  reflected 
back  into  the  jet  plume.  The  reflection  process  is  repeated  many  times  downstream  until  the 
shock/expansion  fan  is  dissipated  by  turbulence.  It  is  these  repeated  reflections  of  the 
shock/expansion  fan  by  the  mixing  layer  of  the  jet  which  give  rise  to  the  quasi-periodic  shock 
cells.  From  this  point  of  view,  the  shock  cell  structure  may  be  regarded  as  disturbances  trapped 
inside  the  jet  by  the  mixing  layer  surrounding  the  jet  column.  In  other  words,  one  may  consider 
the  jet  flow  as  behaving  like  a  waveguide  for  the  disturbances  which  form  the  shock  cell. 

Prandtl  (1904)  appeared  to  be  the  first  to  analyze  the  shock  cell  of  a  supersonic  jet.  He  used  a 
linear  vortex  sheet  jet  model  and  developed  a  partial  solution  to  estimate  the  shock  cell  spacing. 
The  complete  vortex  sheet  shock  cell  solution  for  slightly  imperfectly  expanded  jets  was  given 
later  by  Pack  (1950).  This  solution  has  a  simple  form  in  cylindrical  coordinates  (r,^zi,x).  Let  ps  be 
the  pressure  disturbance  associated  with  the  shock  cells,  then 

00 

A  =  Z  4®,0cos(/r.x)  (1 .23) 

i=\ 

where  H,-  =  2ls.pl (Ji,  0/(r)=Jo(2o^r)/Ji(o^) 


z  =1,2,3,K  . 


Here,  Isp  is  the  static  pressure  difference  inside  and  outside  the  jet  at  the  nozzle  exit.  Jo  and  Ji  are 
the  Bessel  functions  of  order  0  and  1.  Dj  is  the  fully  expanded  jet  diameter  and  A/,  is  the  jet  Mach 
number.  For  a  supersonic  jet  issuing  from  a  convergent-divergent  nozzle  of  design  Mach  number 
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Md{Md=  1.0  for  a  convergent  nozzle)  at  fully  expanded  Maeh  number  A/,,  the  fully  expanded  jet 
diameter  Dj  is  related  to  the  nozzle  exit  diameter  D  by, 
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where  y  is  the  ratio  of  speeifie  heats.  Equation  (1.23)  may  be  eonsidered  as  a  deeomposition  of 
the  shoek  eell  strueture  into  the  waveguide  modes  of  the  jet  flow.  Here  (h/  and  /q  are  the 
amplitude,  eigenfunetion  or  mode  shape  and  wavenumber  of  the  waveguide  mode.  The  shoek 
eell  spaeing,  Ls,  is  given  by  the  longest  wavelength  {n  =  1);  i.e.. 


Itt 


(1.24) 


Obviously,  the  vortex  sheet  shoek  eell  model  is  not  valid  exeept  near  the  nozzle  exit,  where 
the  mixing  layer  of  the  jet  is  thin.  Tam  et  al.  (1985)  extended  the  linear  shoek  eell  solution  to  jets 
with  realistie  mean  flow  profile  using  the  method  of  multiple-seales  expansion.  Figure  1.28 
shows  a  eomparison  between  their  ealeulated  shoek  eell  pressure  distribution  and  the 
measurements  of  Norum  &  Seiner  (1982a).  On  eonsidering  that  the  model  is  linear,  the 
agreement  between  ealeulation  and  measurement  is  remarkably  good  espeeially  in  terms  of  the 
gross  features  sueh  as  the  shoek  eell  spaeing  and  amplitude. 


Figure  1.28.  Comparison  between  ealeulated  and  measured  shoek  eell  axial  pressure 
distribution  at  r/Z)  =  0.38,  Mj=\.%2,  Mrf  =  2.0.  Dark  line  =  ealeulation,  lighter 
line  =  measurement. 
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Dash  et  al.  (1975),  Seiner  et  al.  (1975),  Abdol-Hamid  &  Wilmoth  (1989)  and  others  have 
attempted  to  ealeulate  the  shoek  eells  in  highly  imperfeetly  expanded  supersonie  jets 
numerieally.  They  found  that  a  good  turbulenee  model  was  essential  to  a  good  predietion. 
However,  for  strongly  soreeehing  jets,  all  methods  failed  to  eapture  the  rapid  deeay  of  the  shoek 
eells. 

Tam  (1990)  examined  the  extensive  shoek  eell  measurements  of  Norum  &  Seiner  (1982a) 
and  found  that  for  moderately  imperfeetly  expanded  supersonie  jets  the  first  shoek  was  very 
strong.  However,  downstream  of  the  strong  first  shook,  the  shook  eell  struoture  resembled  that  of 
a  weakly  imperfeetly  expanded  jet.  Based  on  this  observation,  he  suggested  the  use  of  the  linear 
shook  eell  solution  to  model  even  the  shook  eells  of  moderately  imperfeetly  expanded  jets  (not 
valid  for  the  first  shook,  whioh  is  unimportant  as  far  as  broadband  shook  assooiated  noise  is 
oonoerned).  He  developed  a  semi-empirioal  formula  by  whioh  the  initial  amplitude  of  the  linear 
shook  eells  oould  be  estimated. 

1.5.2  Noise  Generation  Mechanism 

Broadband  shook  assooiated  noise  is  generated  by  the  weak  interaotion  between  the 
downstream  propagating  large  turbulenee  struotures  of  the  jet  flow  and  the  quasi-periodio  shook 
eell  struoture  during  the  passage  of  the  former  through  the  latter.  It  is  interesting  to  see  how  suoh 
interaotion  oould  lead  to  the  radiation  of  sound  waves  in  the  upstream  direotions  and  possess 
oharaoteristios  as  desoribed  earlier.  For  simplioity,  a  one-dimensional  model  of  the  large 
turbulenee  struotures-shook  eell  interaotion  is  used  to  explain  the  noise  generation  prooess.  A 
oomplete  three-dimensional  analysis  oan  be  found  in  the  work  of  Tam  (1987). 

Let  Ut  be  the  velooity  fluotuation  assooiated  with  the  large  turbulenee  struotures/instability 
waves.  Within  the  one-dimensional  model  it  may  be  represented  mathematioally  by  a  traveling 
wave  of  the  form 


=  Re[He  (1.25) 

where  Re[  ]  is  the  real  part  of  [  ].  k  and  co  are  the  wavenumber  and  frequenoy  of  the  wave 
disturbanoe.  A  is  the  amplitude.  If  Uc  is  the  oonveotion  velooity  of  the  large  turbulenee 
struotures/instability  waves,  then  k=  cduc-  Now,  in  one  dimension,  the  velooity  fluotuation,  Us 
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associated  with  the  waveguide  mode  of  the  shoek  eells  may  be  represented  by  a  spatially 
periodie  funetion.  Without  loss  of  generality,  let 


=  =  +  (1.26) 

where  k,  is  the  f'  shoek  eell  wavenumber  and  B  is  the  amplitude.  The  eorresponding  shoek  eell 
spaeing  is  equal  to  2p/ki. 

The  interaetion  between  the  downstream  propagating  large  turbulenee  struetures/instability 
waves  and  the  shoek  eells  will  give  rise  to  disturbanees  eonsisting  of  produet  terms  of  equations 
(1.25)  and  (1.26).  The  interaetion  term  most  relevant  to  broadband  shook  assooiated  noise  is 
given  by  the  produet  of  Ut  and  the  seoond  term  on  the  right  side  of  (1 .26).  It  is 

Rej^^  expji^A:  -  -  cot^ .  ( 1 .27) 

(1.27)  represents  a  traveling  wave  with  wave  number  {k-k^.  If  ki  is  slightly  larger  than  k,  then  the 
phase  speed  C  =  a>l{k-ki)  of  the  interaetion  wave  is  negative  and  supersonio.  Aooording  to  the 
wavy  wall  analogy  (see  Figure  1.20)  it  would  immediately  lead  to  Maoh  wave  radiation  in  the 
upstream  direotion. 

The  direotion  of  radiation  q  is  related  to  the  phase  speed  C  and  the  ambient  sound  speed  aoo 
by  the  Maoh  wave  relation;  i.e., 


^00 

C  CO 


=  oos^. 


(1.28) 


Now,  in  terms  of  shook  eell  spaeing  Li  =  Ipiki,  turbulenee  oonveotion  velooity  Uc  and  frequenoy 
/=  coUp,  (1.28)  may  be  rewritten  as 


statement  A:  Approved  for  public 
release;  distribution  is  unlimited. 


50 


AEDC-TR-06-3 


^  rf-,  WC0S6'V  (1-29) 

Ai  1  -  ^ - 1 

\  ^00  J 

As  noted  by  Tam  &  Tanna  (1982),  (1.29)  gives  the  relationship  between  the  spectral  peak 
frequencies  and  the  direction  of  radiation  of  broadband  shock  associated  noise.  This 
characteristic  property  was  discussed  before  and  can  be  seen  in  the  data  of  Figure  1.18.  Since  the 
shock  cells  can  be  considered  as  composed  of  a  number  of  waveguide  modes,  (1.29)  implies  that 
the  total  noise  spectrum  is  made  up  of  a  superposition  of  many  spectral  peaks.  The  spectral  peak 
frequencies  can  be  calculated  by  (1.29). 

The  above  one-dimensional  model  illustrated  the  crucial  noise  generation  processes  of 
broadband  shock  associated  noise.  Basically  this  noise  component  is  generated  by  the 
constructive  scattering  of  the  large  turbulence  structures  of  the  jet  flow  by  the  stationary  quasi- 
periodic  shock  cells  in  the  jet  plume.  Since  the  shock  cells  can  be  conceived  as  made  up  of  a 
superposition  of  waveguide  modes  of  different  wavelengths,  each  mode  scatters  off  sound  in  a 
preferred  direction.  This  results  in  a  multi-peak  noise  spectrum  as  well  as  the  observed 
directional  dependence  of  the  spectral  peak  frequency. 

1.5.3  Stochastic  Model  Theory 

In  a  series  of  papers  (1987,  1990,  1991b,  1992),  Tam  developed  a  semi-empirical  stochastic 
model  theory  for  the  prediction  of  broadband  shock  associated  noise.  In  formulating  the  theory, 
Tam  divided  the  flow  variables  into  four  components,  namely,  the  mean  flow,  the  large 
turbulence  structures/instability  waves,  the  shock  cells  and  the  disturbances  arising  from  the 
interaction  between  the  large  turbulence  structures/instability  waves  and  the  shock  cells.  The  last 
of  the  four  components  contained  the  broadband  shock  associated  noise.  Tam  showed  that  this 
component  could  be  found  by  solving  a  nonhomogeneous  boundary  value  problem.  The 
nonhomogeneous  terms  were  the  product  terms  of  the  large  turbulence  structures/instability 
waves  and  the  shock  cells.  Tam  (1987)  solved  this  problem  formally.  However,  to  make  use  of 
the  formal  solution  for  noise  prediction,  extensive  numerical  computations  were  required.  To 
avoid  such  lengthy  computations,  Tam  introduced  a  similarity  spectrum  to  approximate  the 
disturbances  generated  by  the  turbulence  shock  cell  interaction.  The  similarity  spectrum  was 
based  on  the  assumption  that  the  noise  source  spectrum  (the  large  turbulence  structures)  had  no 
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intrinsic  length  and  time  seales.  When  this  is  true  the  speetrum  must  be  a  funetion  of fxlUj  where 
/is  the  frequeney,  x  is  the  downstream  distanee  and  Uj  is  the  fully  expanded  jet  veloeity.  The 
semi-empirieal  theory  of  Tam  eontains  four  parameters,  two  of  whieh  ean,  in  prineiple,  be 
ealeulated  by  means  of  the  instability  wave  theory.  The  unknown  parameters  were  determined 
empirieally  by  Tam  (1987,  1990). 

Let  Uc,  Uj  and  w/be  the  eonveetion  veloeity  of  the  instability  waves,  the  jet  veloeity,  and  the 
flight  veloeity,  respeetively.  If  /  is  the  total  temperature  of  the  jet  and  Too  is  the  ambient 
temperature,  then  Uc  may  be  estimated  by  the  following  formula: 


M,. 


fT  Y 

0.7-0.06|^^-1.0j 


Uf 

0.7  +  0.3^, 

M, 


1 - ^ 


Uf  T 

H - ,  >  1 

V  ^Jj  Uj 


T 

^<1 


(1.30) 


The  half-width,  L,  of  the  similarity  noise  souree  may  be  estimated  by  the  empirieal  formula. 


3.3 


L  = 


Vs.oj 


f 

1.0  + 

1.114-0.367/ 

\Mj. 

D: 


(1.31) 


where  Xc,  the  core  length  of  the  jet,  ean  be  ealeulated  by 


A3  +  \.2M], 

T,>T^ 

< 

r  f 

4.3  +  1.2M"+  1 

V 

yY 

Tj 

,  T,<T^ 

(1.32) 


where  7)  is  the  temperature  of  the  jet,  and  Dj  is  the  fully  expanded  jet  diameter  related  to  the 
nozzle  exit  diameter  D  by 
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1  + 

(/-I) 

2 

Mj 

D 

1  + 

M] 

2 

fr+i) 

4(7-1) 


kMjJ 


(1.33) 


where  Md  is  the  nozzle  design  Maeh  number.  The  shoek  eell  spaeing,  Lm,  of  hot  jets  is  related  to 
wavenumber  km  by  Lm  =  27rlkm-  The  wave  number  of  the  shoek  eell  waveguide  mode  km  may 
ealeulated  by, 


„(hot  jet  atM/)= 


^old  jet  atM^  =  0 

) 

1.0  + 

^  0.2547,3 

U.O  1.Z,  rp  J 

V  ^00  J 

(1.34) 


with  the  parameters  of  the  shoek  eell  strueture  determined.  A  broadband  shook  assooiated  noise 
power  speotrum  formula  for  hot  supersonio  jets  in  flight  in  the  nozzle  fixed  ooordinates  is 
derived  by  Tam  (1995)  as, 


5 


cL^A^A^platM] 


R 


y/y 


fD, 


1  + 


(r-1) 


I 


m=l 


■exri 


/  A 

2 

Jm  1 

1/  J 

V 

l  +  M. 


\Mi 


M j-  ^  -  M j  sin^  +  cos 


^  -  M^j-  sin^  y?^ 


21n2 


(1.35) 


where  c  =  2.65  x  IO^^q  (So  =  0.35),  Aj  =  nD^.  / 4 ,  and  is  the  ambient  speed  of  sound.  J\  is  the 

Bessel  funotion  of  order  one,  and  cr„,  is  the  tn^  zero  of  the  Bessel  funotion  of  order  zero.  In  (1.35) 
the  faotor  fDj/uj  in  the  denominator  is  to  be  replaoed  by  5*0  whenever  it  is  less  than  5*0  and 


f  = 

J  m 


u  k 

cm 


l  +  M, 


sin^  1//J  +  oos^ 
^  -  M^  sin^  yj 


> 


(1.36) 
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In  formulating  the  broadband  shock  associated  noise  theory  it  has  always  been  assumed  that 
the  convergent-divergent  nozzle  of  is  perfectly  manufactured.  That  is,  when  operating  at  a 
fully  expanded  Mach  number  Mj  =  Md,  there  is  no  shock  cell  structure  in  the  jet  plume.  In 
practice  this  is  not  always  true,  especially  for  hot  jets.  Because  of  various  nozzle  imperfections,  a 
residual  shock  cell  structure  invariably  exists  so  that  even  when  the  jet  is  operating  nominally  at 
the  design  Mach  number,  there  is  still  broadband  shock  associated  noise.  To  account  for  this 
deviation  from  an  ideal  condition,  one  may  incorporate  empirically  a  residue  shock  cell  strength 
8  in  the  parameter  of  (1.35)  as  follows: 


A^  =  { 


^  / 

- 

-  /  1  +  3 

K-M;)  +  £| 


1  +  ^M' 


(underexpanded  j  ets) 


1  +  6 


I 


(1.37) 


(overexpanded  j  ets) 


where  ^  =  0  corresponds  to  a  C-D  nozzle  without  imperfections.  Based  on  a  number  of  numerical 
test  cases,  it  is  found  that  a  value  of  ^  =  0.03  gives  generally  satisfactory  results  compared  with 
experimental  measurements.  It  should  be  noted  that  the  residual  shock  cell  strength  is  important 
only  when  the  jet  is  operating  closed  to  its  nominal  condition. 

Extensive  comparisons  between  the  semi-empirical  noise  prediction  formula  and 
experimental  measurements  have  been  carried  out.  For  all  the  cases  compared,  agreement  good 
enough  for  engineering  application  was  found.  This  was  true  regardless  of  whether  the  jet  was 
hot  or  cold  and  whether  the  jet  was  underexpanded  or  overexpanded.  In  these  tests,  the  nozzle 
design  Mach  number  and  the  fully  expanded  jet  Mach  number  covered  the  range  of  1.0  to 
slightly  over  2.0.  Figure  1.29  shows  typical  calculated  noise  spectra  and  the  corresponding 
experimental  measurements.  As  can  be  seen,  except  for  the  dips  at  small  inlet  angle  directions, 
there  is  good  agreement  between  the  calculated  and  the  measured  spectra. 
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The  semi-empirical  theory  is  also  capable  of  predicting  the  near  field  broadband  shock 
associated  noise.  This  is  important  from  the  standpoint  of  sonic  fatigue  since  structures  near  the 
tail  part  of  a  jet  aircraft  are  exposed  to  constant  broadband  shock  associated  noise  radiation. 
Calculated  near  field  noise  contours  and  their  comparisons  with  the  measurements  of  Yu  (1970) 
are  provided  in  the  work  of  Tam  (1987).  The  contour  shapes  and  distributions  are  shown  to  agree 
well  with  measurements. 


Figure  1.29.  Comparisons  between  calculated  broadband  shock  noise  spectra  and 
measurements  of  Norum  &  Seiner  (1982a).  Mj  =  1.67,  Md  =1.0 
The  original  effort  in  developing  semi-empirical  broadband  shock  associated  noise  prediction 

method  was  primarily  for  circular  jets.  Later  Tam  and  Reddy  (1996)  extended  the  approach  to 

rectangular  jets.  The  change  in  geometry  requires  modification  to  the  shock  cell  model  as  well  as 

to  the  large  turbulence  structures/instability  wave  model.  An  extra  parameter,  namely,  the  aspect 
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ratio  of  the  jet,  has  to  be  introdueed.  The  predietions  of  the  extended  theory  were  shown  to  agree 
fairly  well  with  measurements.  However,  the  data  base  is  limited.  Further  testing  of  the  semi- 
empirical  theory  would  provide  added  confidence  in  its  accuracy. 


1,6  Screech  Tones 

Screech  tones  are  very  sensitive  to  changes  in  the  surrounding  environment.  This  makes  them 
somewhat  unpredictable.  For  instance.  Norum  (1983)  found  that  the  screech  intensity  increased 
by  lOdB  when  a  thin  lip  nozzle  was  replaced  by  one  with  a  thicker  lip.  Moreover,  it  has  been 
reported  by  Seiner  et  al.  (1986)  that  the  flapping  plane  of  a  screeching  jet  was  found,  for 
unknown  reasons,  to  have  rotated  by  almost  90  degrees  when  the  same  experiment  was  repeated 
a  month  later  in  the  same  facility  by  the  same  investigators. 

1,6,1  Tone  Generation  Mechanism 

Powell  (1953a,b)  was  the  first  to  report  the  observation  of  screech  tones.  He  correctly 
suggested  that  it  was  an  acoustic  feedback  phenomenon.  Over  the  years,  our  understanding  of  the 
details  of  the  screech  feedback  loop  has  undergone  significant  changes.  Figure  1.30  shows  the 
present  day  view  of  the  various  components  of  the  feedback  loop.  Near  the  nozzle  lip  where  the 
jet  mixing  layer  is  thin  and  most  receptive  to  external  excitation,  acoustic  disturbances  impinging 
on  this  area  excite  the  intrinsic  instability  waves  of  the  jet  flow.  The  amplitude  of  the  excited 
instability  is  small  near  the  nozzle  exit.  However,  as  the  instability  wave  propagates  downstream 
it  extracts  energy  from  the  mean  flow  and  grows  rapidly  in  amplitude.  After  propagating  a 
distance  of  four  to  five  shock  cells,  the  instability  wave,  having  acquired  a  large  enough 
amplitude,  interacts  with  the  quasi-periodic  shock  cells  in  the  jet  plume.  The  unsteady  interaction 
generates  acoustic  radiation.  For  the  same  reason  as  in  the  case  of  broadband  shock  associated 
noise,  the  acoustic  radiation  is  primarily  in  the  upstream  direction  as  illustrated  in  Figure  1.30. 
The  feedback  acoustic  waves  propagate  upstream  outside  the  jet.  Upon  reaching  the  nozzle  lip 
region,  they  excite  the  shear  layer  of  the  jet  leading  to  the  generation  of  new  instability  waves.  In 
this  way,  the  feedback  loop  is  closed. 
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Figure  1.30.  Schematic  diagram  of  the  screech  tone  feedback  loop. 

1.6.2  Fundamental  Screech  Tone  Frequency 

Tam  et  al.  (1986)  suggested  that  the  weakest  link  of  the  screech  feedback  loop  is  at  the 
nozzle  lip  where  feedback  acoustic  waves  excite  the  instability  wave.  Their  contention  is  that 
unless  the  acoustic  waves  are  sufficiently  strong,  the  excitation  process  will  not  be  able  to 
generate  instability  waves  at  a  large  enough  amplitude  (above  the  background  noise)  to  maintain 
the  feedback  loop.  To  maintain  the  feedback  process,  it  is,  therefore,  essential  that  the  direction 
of  maximum  acoustic  radiation  generated  by  the  passage  of  the  instability  wave  through  the 
shock  cells  of  the  jet  be  in  the  direction  of  the  nozzle  lip.  The  relation  between  the  sound 
frequency  and  the  direction  of  radiation  is  given  by  (1.29).  By  setting  6=  n,  we  achieve  a 
formula  by  which  the  fundamental  screech  frequency  can  be  estimated.  By  adapting  the 
empirical  results  of  Uc  =  0.76)  (6)  =  fully  expanded  jet  velocity)  and  the  shock  cell  spacing  at 
four  to  five  shock  cells  downstream  to  be  about  20%  smaller  than  that  given  by  the  vortex  sheet 
shock  cell  model  (1.24),  Tam  et  al.  (1985)  derived  the  following  screech  tone  frequency  formula 
for  a  jet  of  Mach  number  A/;. 
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where  fs  is  the  fundamental  sereeeh  frequeney.  Tr  is  the  total  temperature  of  the  jet  and  is  the 
ambient  temperature.  Dj  is  the  fully  expanded  jet  diameter  (see  equation  1.33). 


) 


1  +  0.1  M. 


f 


(1.38) 


Figure  1.31.  Comparisons  between  measured  and  ealeulated  sereeeh  tone 
frequencies  at  different  total  temperature  J). 
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Figure  1.31  shows  a  comparison  between  (1.38)  and  the  measurements  of  Rosfjord  &  Toms 
(1975)  at  three  jet  temperatures.  In  the  experiment  a  convergent  nozzle  {Md  =  1.0)  was  used.  The 
agreement  between  prediction  and  measurements  is  quite  good.  It  is  better  at  high  temperature 
and  high  pressure  ratio. 
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Figure  1.32.  Rectangular  jet  with  aspect  ratio  hj/bj. 

Tam  (1988)  used  similar  ideas  to  derive  a  screech  frequency  formula  for  rectangular  jets. 
Consider  a  vortex  sheet  bounded  rectangular  jet  with  fully  expanded  width,  bj,  and  height,  hj,  as 
shown  in  Figure  1.32.  For  weak  shocks,  the  governing  equations  are  the  linearized  Euler 
equations. 


dp 


ds  ^ 


p  =  a  p. 


(1.39(a)) 


(1.39(b)) 


(1.39(c)) 


By  elimination,  a  single  governing  equation  for  p  is 

V^p-M]^  =  0  (1.40) 

^  dK 

where  Mj  is  the  Mach  number.  At  the  vortex  sheet,  linearized  form  of  the  dynamic  boundary 
condition  (pressure  balance)  yields 


p  =  0  at  the  vortex  sheet. 


(1.41) 
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It  will  be  assumed  that  the  pressure  at  the  nozzle  exit,  pe,  is  uniform  aeross  the  entire  eross 
seetion  of  the  jet  and  differs  from  the  ambient  pressure,  pa,  by  an  amount  ISp.  It  will  further  be 
assumed  that  at  the  nozzle  exit  the  flow  veloeity  is  paralleled  to  the  x-axis  or  the  direetion  of  jet 
flow.  Thus  the  initial  eonditions  for  the  shoek  eell  solution  at  the  nozzle  exit  x  =  0  are, 

p  =  bq)  and  =  0 .  (1-42) 

The  solution  of  (1.40)  satisfying  boundary  eondition  (1.41)  and  nozzle  exit  eondition  (1.42) 
is 


p{x,y,z)=Y,'Y. 

n=\  m=\ 


n  Yi 

— ^1  -  cosnTTjy  -  eos 
nmn 


X  .  HTiy  .  mTiy  , 

— cos  k^jc 


(1.43) 


where 


n,m  =  1,2,3,K  . 


(1.44) 


The  smallest  eigenvalue  or  wave  number  is  ku.  Thus  the  shoek  eell  spaeing  Ls  is  equal  to 


lihjibj  <  1/4,  then  (1.45)  may  be  approximated  by 


4=2(m;-i)%, 


(1.45) 


(1.46) 


whieh  is  the  shoek  eell  spaeing  of  a  two-dimensional  jet  with  thickness  hj. 
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To  compare  (1.46)  with  experiments,  it  is  neeessary  to  relate  the  fully  expanded  jet  width  bj 
and  height  hj  to  the  nozzle  exit  dimensions.  Tam  and  Tanna  (1982)  pointed  out  that  the  effeetive 
size  of  a  fully  expanded  supersonic  jet  is  not  the  same  as  that  at  the  nozzle  exit.  The  differenee  is 
a  funetion  of  jet  operating  eonditions  and  eould  sometimes  be  signifieant.  To  aeeount  for  the  size 
differenee  let  and  Aj  be  the  eross-seetional  area  of  a  jet  at  the  nozzle  exit  and  that  of  the  fully 
expanded  jet,  Md  and  Mj  be  the  nozzle  design  Maeh  number  and  the  fully  expanded  Maeh 
number.  From  one-dimensional  gas  dynamies,  it  is  easy  to  see  that  Aj  and  Ad  are  related  by  (see 
the  derivation  provided  by  Tam  (2005)), 


Ms 

m' 


(1.47) 


Let  b  and  h  be  the  nozzle  exit  dimensions;  then,  by  assuming  the  differenees  between  the  width 
and  height  of  the  fully  expanded  jet  and  b  and  h  are  small,  it  is  simple  to  find. 


h 

h 


(1.48(a)) 


h  +  b 


+  1. 


(1.48(b)) 


In  the  ease  b  »  h,  the  shoek  eell  spaeing  as  given  by  (1.46),  (1.47)  and  (1.48)  reduees  to. 
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Figure  1.33.  Comparison  between  measured  and  ealeulated  shoek  eell  spaeings  of 
reetangular  jets.  #,  Powell  (1953),  b/h  =  5.S3;  O,  Hammitt  (1961),  blh=  12.75; 
- ,(1.46) 

Figure  1.33  shows  eomparisons  between  (1.49)  and  the  measured  data  of  Powell  (1953c)  and 
Hammit  (1961).  There  is  good  agreement  over  a  wide  range  of  jet  Mach  number. 

According  to  the  weakest  line  screech  feedback  tone  model  of  Tam,  Seiner  and  Yu  (1986), 
the  screech  frequency  may  be  obtained  by  setting  0  =  tt  and  Ls  =  Ln  in  (1.29).  This  gives  the 
following  screech  tone  frequency  formula  for  rectangular  jets. 


U: 


+  1 


V  aj 


(mJ-i) 


r+1 

2(7-1) 


M,. 


-1 


b 

^b  +  h^^ 


(1.50) 


For  cold  choked  jets  {Md  =  1.0)  with  large  aspect  ratio  (1.50)  simplifies  to. 
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^1+  2 

(1.51) 


ujuj  =  0.7  has  been  used  to  derive  (1.51). 

Figure  1.34  shows  a  eomparison  of  the  numerieal  results  of  (1.51)  and  the  experimental 
measurements  of  Raman  (1996).  There  is  exeellent  agreement  between  predietions  and 
measurements  over  the  Maeh  number  range  of  1.1  <M,  <  1.8. 


Figurel.34.  Variation  of  the  sereeeh  frequeney  with  fully  expanded  jet  Maeh 
number  for  a  supersonie  jet  from  a  eonvergent  reetangular  nozzle;  aspeet 
ratio  =  5,  regular  exit  geometry;  O,  measurements  and - ,  theory. 
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Some  military  jets  use  beveled  rectangular  nozzles  as  shown  in  Figure  1.35.  The  shock  cell 
structure  of  jets  issuing  from  this  type  of  nozzle  is  highly  complex.  However,  beyond  the  most 
downstream  edge  of  the  nozzle,  the  shock  cell  structure  is  quasi-periodic.  This  is  the  region 
where  the  feedback  acoustic  waves  of  the  screech  feedback  loop  are  generated.  The  physical 
situation  is  very  similar  to  that  of  a  jet  with  regular  straight  nozzle  exit.  The  major  difference  lies 
in  that  the  dominant  wave  guide  modes  that  make  up  the  quasi-periodic  shock  cell  structure  are 
not  necessarily  the  lowest-order  modes.  Moreover,  it  is  possible  that  more  than  one  mode  is 
dominant. 


Figure  1.35.  Spark  Schlieren  image  of  the  shock  cell  structure  of  a  Mach  1.4, 
aspect  ratio  5,  rectangular  jet  from  a  single-beveled  convergent  nozzle. 

As  far  as  screech  tone  feedback  loop  is  concerned,  the  weakest  link  theory  is  still  applicable. 
However,  one  expects  the  resulting  screech  tones  to  be  more  complex.  This  is  because  they  are 
related  to  the  higher-order  modes;  also,  more  than  one  feedback  loop  may  be  operating  at  the 
same  time.  On  following  the  weakest  link  theory,  Tam,  Shen  and  Raman  (1997)  proposed  that 
the  screech  frequencies  can  be  calculated  by  replacing  Li  in  (1.29)  by  Iniknm  where  knm  are  given 
by  (1.44).  For  large  aspect  ratio  jets  {bj  »  kj)  only  the  m  =  \  modes  are  relevant.  On  restricting 
consideration  to  m  =  1  modes,  the  following  screech  frequency  formula  is  obtained; 
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(1.52) 


where  n  =  1,2,3,...  . 

The  sereeeh  frequeneies  of  supersonie  jets  from  beveled  reetangular  nozzles  have  been 
measured  by  Raman  (1996).  The  frequeney  as  a  funetion  of  jet  Maeh  number  has  a  fairly 
eomplex  pattern.  Figure  1.36  shows  the  measured  sereeeh  tone  frequeneies  as  a  funetion  of  the 
jet  Maeh  number  for  a  supersonie  jet  from  a  single-beveled  eonvergent  reetangular  nozzle.  The 
data  do  not  seem  to  form  a  eontinuous  curve.  There  are  frequency  jumps  at  several  jet  Mach 
numbers.  Also  shown  are  the  screech  frequencies  of  (1.52)  with  spanwise  mode  number  n=  1,2, 
3,  and  4.  As  can  be  seen,  the  data  correlate  well  with  the  theoretical  frequency  curves  over  the 
entire  range  of  measured  Mach  numbers.  Figure  1.37  shows  a  comparison  of  the  measured 
screech  frequencies  and  the  frequencies  of  (1.52)  for  rectangular  jets  from  a  double-beveled 
convergent  nozzle.  The  screech  frequencies  agree  well  with  the  n  =  5  and  6  frequency  curves. 

Figure  1.38  shows  the  variation  of  the  screech  frequency  with  Mach  number  for  the 
supersonic  jet  from  a  C-D  rectangular  nozzle  with  a  regular  exit  geometry.  Based  on  his 
spanwise  phase  data,  Raman  was  able  to  identify  three  distinct  screech  modes.  Their  frequencies 
fit  very  well  with  the  frequencies  of  the  m  =  1,  3,  and  5  modes  of  (1.52).  Figure  1.39  provides  the 
measured  screech  frequency  jet  Mach  number  relation  of  the  supersonic  jet  from  a  single- 
bevelled  C-D  rectangular  nozzle.  According  to  the  phase  data  there  are  three  modes;  i.e.,  modes 
1,  lllA,  and  lllB.  Also  shown  are  the  screech  frequency  curves  for  m  =  1,  2,  3,  and  4  modes 
according  to  (1.52).  The  data  appear  to  agree  well  with  the  frequency  curves  of  the  m  =  1,  2,  and 
4  modes.  Finally,  Figure  1.40  shows  the  measured  screech  frequency  variation  with  jet  Mach 
number  of  a  supersonic  jet  from  a  C-D  double-beveled  rectangular  nozzle.  Two  modes  have  been 
identified.  Their  frequencies  are  in  reasonable  agreement  with  the  n  =  6  and  7  frequency  curves 
of(1.52). 

Equation  (1.52)  predicts  many  screech  frequency  bands.  However,  not  all  of  the  tones  of  the 
frequency  bands  are  excited.  It  is  believed  that  the  selection  of  the  excited  frequencies  is 
determined  largely  by  two  factors.  These  are  the  spatial  growth  rates  of  the  instability  waves  of 
the  jet  flow  and  the  strength  of  the  waveguide  modes  of  the  shock  cell  structure. 
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Figure  1.36.  Variation  of  the  sereeeh  frequeneies  with  fully  expanded  jet  Mach 
number  for  a  supersonic  jet  from  a  convergent  rectangular  nozzle  with  aspect 

ratio  =  5,  single-beveled  exit  geometry:  A,  measurements.  (1.52): - ,  n= 

- ,n  =  2;  —  •  — ,  n  =  3;  and - ,  n  =  4. 
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Figure  1.37.  Variation  of  the  screeeh  frequeneies  with  fully  expanded  jet  Maeh 
number  for  a  supersonic  jet  from  a  convergent  rectangular  nozzle  with  aspect 

ratio  =  5,  double-beveled  exit  geometry:  □,  measurements.  (1.52): - ,  n= 

- ,n  =  2;  —  •  — ,  n  =  3; - ,  n  =  4;  —  •  •  — ,  n  =  5;  and . ,n  =  6. 
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Figure  1.38.  Variation  of  the  sereeeh  frequencies  with  fully  expanded  jet  Mach 
number  for  a  supersonic  jet  from  a  convergent-divergent  rectangular  nozzle  with 
aspect  ratio  =  5,  Mj  =  1.4,  and  regular  exit  geometry.  Measurements:  O,  mode  11; 
A,  mode  12;  and  □,  mode  11.  (1.52): - ,  n  =  1; - ,n  =  3;  —  •  — ,  n  =  5. 
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Figure  1.39.  Variation  of  the  sereech  frequencies  with  fully  expanded  jet  Mach 
number  for  a  supersonic  jet  from  a  convergent-divergent  rectangular  nozzle  with 
aspect  ratio  =  5,  Md=  1.4,  and  single-beveled  exit  geometry.  Measurements:  O, 

mode  1;  A,  mode  llIA;  and  □,  mode  lllB.  (1.52): - ,  n  =  1; - ,  n  =  2; 

—  •  — ,n  =  3;  and - ,  n  =  4. 
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Figure  1.40.  Variation  of  the  sereeeh  frequeneies  with  fully  expanded  jet  Mach 
number  for  a  supersonic  jet  from  a  convergent-divergent  rectangular  nozzle  with 
aspect  ratio  =  5  and  double-beveled  exit  geometry.  Measurements:  O,  mode  1  and 
A,  mode  11.  (1.52): - ,n  =  6  and . ,n  =  l. 
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Chapter  2.  The  Acoustic  Environment  Formed  by  a  Test  Cell 

An  engine  test  eell,  invariably,  provides  a  eonfined  environment  for  aeoustie  waves.  The 
eonstant  refleetions  of  acoustic  waves  by  the  solid  walls  and  opened  ends  of  an  engine  test  cell 
lead  to  the  formation  of  standing  waves.  These  standing  waves,  often  referred  to  as  normal 
modes,  allow  the  accumulation  of  wave  energy.  This  results  in  large  amplitude  oscillations  and 
resonances.  For  test  cells  with  relatively  simple  geometry,  the  resonance  frequencies  can  be 
found  analytically.  If  the  geometry  is  significantly  complicated,  CAA  methods  may  have  to  be 
used. 

2,1  Ejector  Type  Engine  Test  Cells 

Ground  level  test  cells  are  effectively  large  ejectors  with  open  ends.  The  geometry  may  be 
circular  or  rectangular  as  shown  in  Figure  2.1. 


X 


L 


(a)  Cylindrical  geometry 


Figure  2.1.  Ejector  type  ground  level  engine  test  cells. 
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Figure  2.1.  Ejector  type  ground  level  engine  test  cells. 


Let  us  assume  that  there  is  a  uniform  mean  flow,  mq,  in  the  axial  or  the  x-direction.  The 
governing  equations  for  the  acoustic  disturbances  are  the  linearized  Euler  equations. 


(d^  d^'\ 


-V/7 


(2.1) 


dp  dp 

■^  +  “o  ^  +  7P^  •  V  =  0 


di 


(2.2) 


where  /jo,  />o  are  the  density  and  pressure  of  the  mean  flow  and  yis  the  ratio  of  specific  heats  of 
the  gas.  Upon  eliminating  V,  the  governing  equation  for  p  is. 


(d 


dK.) 


p  -  alV^p  =  0 


(2.3) 


1/2 

where  ao  =  is  the  speed  of  sound. 

2.1.1  Cylindrical  Test  Cells 
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For  cylindrical  test  cells,  the  natural  coordinates  to  use  are  the  eylindrieal  eoordinates  (r,(j)^ 
as  shown  in  Figure  2.1.  Equation  (2.3)  beeomes, 


d'^  2 

p-a^ 


^p  1 


+  - 


■  + 


^ p  p 


r  dr  d  d(/P'  ) 


=  0. 


(2.4) 


The  normal  mode  solutions  have  the  form, 

=  (2.5) 


where  Re{  }  is  the  real  part  of  k,  the  axial  wave  number  and  m,  an  integer,  is  the  azimuthal 
mode  number;  co  is  the  angular  frequeney.  Substitution  of  (2.5)  into  (2.4)  gives. 


p  1  dp 
dd  ^  r  dr  ^ 


(2.6) 


The  solution  of  (2.6),  whieh  is  bounded  at  r  =  0  is. 
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where  Jm{  )  is  the  m*  order  Bessel  funetion. 

The  boundary  eondition  at  the  eylindrieal  surface  is  Vr  =  0.  By  (2.1),  this  is  equivalent  to. 


r  =  -  ^  =  0 

2’  dr 


Let  Amn  be  the  ri^  zero  of  the  derivative  Jm,  i.e., 

n  =  1,2,3,. 


(2.8) 


(2.9) 
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Thus,  substitution  of  solution  (2.7)  into  boundary  condition  (2.8)  yields 


or 


([-Mi'y- +imq 


CO 

a  ok 


4/l„ 


CO 

~  , 


=  0 


(2.10) 


where  Mo  =  uo/ao  is  the  mean  flow  Mach  number.  Equation  (2.10)  relates  k  and  co.  Such  a 
relationship  is  called  the  dispersion  relation.  Dispersion  relation  governs  wave  propagation  in  a 
duct. 

On  solving  (2.10)  for  k,  one  finds  two  solutions.  They  are, 


(2.11) 


where  m  =  0,1,2,...;  n  =  1,2,3,...  .  By  (2.5),  (2.7)  and  (2.11)  the  full  solution  is. 


/»(r,^i^,x,t)=Rej 


\  D  j 


(2.12) 


Now,  at  the  two  open  ends  of  the  test  eell,  the  perturbation  pressure  is  small.  Thus  the 
appropriate  boundary  conditions  are 

x  =  0,  X  =  L,  p  =  0.  (2.13) 

On  imposing  (2.13)  on  (2.12),  it  is  straightforward  to  find, 

(2.14a) 

and 

=1.  (2.14b) 


statement  A:  Approved  for  public 
release;  distribution  is  unlimited. 


74 


AEDC-TR-06-3 


The  solution  of  (2.14b)  is 

By  (2.1 1),  this  condition  when  written  out  in  full  is, 


Upon  solving  for  co,  the  normal  mode  frequencies  are  found  to  be, 


(2.15) 


2.1.2  Rectangular  Test  Cells 

For  test  cells  with  a  rectangular  geometry  of  dimensions  W  by  H  by  L  as  shown  in  Figure 
2.1,  a  similar  analysis  as  above  leads  to  a  dispersive  relation  in  the  form 


(2.16) 


(2.17) 


2,2  High-Altitude  Engine  Test  Cells 

A  high-altitude  engine  test  cell  forms  a  closed  environment  by  itself  The  geometry  of  some 
altitude  test  cells  is  fairly  complex.  The  acoustic  resonance  frequencies  of  these  test  cells  cannot 
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be  written  out  in  simple  analytieal  forms.  They  ean,  however,  be  found  computationally  by  CAA 
methods. 

In  some  high-altitude  test  cells,  a  circular  diffuser  piece  as  shown  in  Figure  2.2  is  sometimes 
used.  Such  a  diffuser  has  resonance  frequencies  of  its  own.  They  can  be  found  as  follows. 

!'■ 


D 

-  L  - >■ 

Front  view  Cross  section 


Figure  2.2.  Schematic  diagram  of  a  diffuser  with  one  open  end  and  one  nearly  closed  end. 

Let  us  approximate  the  diffuser  as  a  short  cylindrical  shell  with  one  open  and  one  closed  end. 
The  solution  inside  the  cylindrical  diffuser  is  still  given  by  (2.12).  But  now,  the  end  conditions 
are. 


X  =L,  p  =  0,  (open  end) 
dp 

X  =  0,  'dx~^’  (closed  end). 

Upon  imposing  boundary  condition  (2.18)  on  (2.12)  it  is  easy  to  find 

=0. 

Similarly,  the  imposition  of  boundary  condition  (2.19)  on  (2.12)  gives 


klA  +  klB  =  0 . 


(2.18) 

(2.19) 


(2.20) 
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By  means  of  this  relation  and  (2.20),  it  is  straightforward  to  find, 


=  k- 

mn  mn 


(2.21) 


With  one  end  closed,  the  effective  mean  flow  velocity  should  be  small.  Thus  it  is  reasonable  to 
set  Mo  =  0  in  calculating  .  This  leads  to 


C  =  = 


®  42? 


(2.22) 


Thus  (2.21)  becomes 
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The  solution  of  (2.23)  is 
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where  i  is  an  integer.  On  solving  for  o),  it  is  found. 


COt  =  Cln 

i,m,n  0 
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(2.24) 


cOf^mn  as  given  by  (2.24)  are  the  normal  mode  frequencies  of  the  cylindrical  diffuser.  They  are  also 
the  acoustic  resonance  frequencies 
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Chapter  3.  Interaction  Noise 

The  insertion  of  a  jet  inside  an  engine  test  eell  may  lead  to  noise  and  tones  that  are  generated 
by  the  coupling  of  the  jet  to  the  acoustic  environment  of  the  test  cell.  Such  interaction  could  lead 
to  resonances  and  super  resonance.  The  energy  source  of  these  resonances  is  primarily  the 
instabilities  of  the  jet  flow.  These  instabilities  are  excited  by  the  acoustic  modes  of  the  test  cell. 
The  excitation  process  is  generally  referred  to  as  receptivity.  Jets  are  unstable  only  over  a  band 
of  frequencies.  The  growth  rate  of  instability  waves  at  different  frequencies  within  the  band  can 
vary  significantly.  When  the  most  amplified  instability  wave  frequency  matches  that  of  a  duct 
mode  of  the  test  cell,  very  high  amplitude  resonance,  called  super  resonance,  could  occur.  This  is 
most  undesirable  and  destructive.  It  may  render  severe  damage  to  the  test  facility. 


3.1  Instabilities  of  High-Speed  Jets 

A  jet  is  intrinsically  unstable.  That  is,  small  perturbations  imposed  on  a  jet  will  have  the 
natural  tendency  to  grow  to  a  substantial  amplitude.  A  comprehensive  study  of  the  instabilities  of 
high-speed  jets  has  been  carried  out  by  Tam  and  Hu  (1989).  They  found  that  these  jets  could 
support  three  families  of  instability  waves.  These  instabilities  have  been  observed  by  Oertel 
(1979,  1980,  1982).  The  three  families  of  instabilities  are  the  Kelvin-Helmholtz  instability,  the 
supersonic  and  subsonic  instability  waves.  The  supersonic  and  subsonic  instability  waves 
propagate  with  phase  velocities  that  are  supersonic  and  subsonic  relative  to  ambient  speed  of 
sound.  Of  importance  to  engine  test  cell  applications  are  the  Kelvin-Helmholtz  instability  waves. 
These  instability  waves  play  a  crucial  role  in  providing  energy  to  drive  the  feedback  loop  of  jet 
screech  tones,  in  the  generation  of  Mach  wave  radiation  and  broadband  shock  noise,  and  in 
providing  an  energy  source  responsible  for  the  occurrence  of  many  engine  test  cell  resonances. 

Before  embarking  on  an  analysis  of  the  instability  wave  modes  of  high-speed  jets,  it  is  found 
useful  first  to  examine  physically  why  these  jets  can  support  three  distinct  families  of  waves. 
Ackeret  (see  Liepmann  &  Puckett,  1947;  Papamoschou  &  Roshko,  1986)  was  the  first  to  provide 
a  physical  explanation  of  the  mechanism  responsible  for  the  Kelvin-Helmholtz  instability  in  a 
thin  shear  layer.  Here  it  will  be  shown  that  Ackeref  s  physical  argument  can  be  extended  to  the 
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case  of  high-speed  jets  to  explain  why  these  jets  can  sustain  three  elasses  of  waves.  For  the  sake 
of  eompleteness,  Ackeret’s  explanation  of  the  Kelvin-Helmholtz  instability  meehanism  will  be 
briefly  reviewed. 

(a) 


u 


(b) 


c 


u  -  c 


Figure  3.1.  Kelvin-Helmholtz  instability  mechanism  a  subsonic  Mach  number. 

+,  high  pressure  region;  low  pressure  region,  (a)  Stationary  frame  of  referenee. 

(b)  Wave  frame  of  referenee. 

Consider  a  two-dimensional  vortex  sheet  separating  a  fluid  at  rest  and  a  fluid  moving  at  a 
subsonie  veloeity  U.  It  will  be  assumed  that  the  vortex  sheet  is  deformed  by  a  Kelvin-Helmholtz 
wave  with  a  phase  veloeity  c  as  shown  in  Figure  3.1(a).  Ackeret  suggested  that  one  should  view 
the  flow,  not  in  the  stationary  frame  of  referenee,  but  in  a  moving  frame  traveling  with  the  phase 
veloeity  c  of  the  wave.  In  this  wave  frame  the  flow  is  as  shown  in  Figure  3.1(b).  Now  for  the 
flow  above  the  vortex  sheet,  the  vortex  sheet  may  be  regarded  as  a  wavy  wall.  So  within  a  quasi¬ 
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steady  approximation,  the  flow  is  that  of  a  uniform  flow  past  a  wavy  wall.  The  solution  of  this 
problem  is  well  known  (see,  e.g.,  Liepmann  &  Roshko,  1957,  Ch.  8).  At  subsonie  Maeh  number 
the  pressure  is  lowest  at  the  erests  of  the  wavy  wall  and  highest  at  the  troughs.  Similar 
eonsideration  may  also  be  applied  to  the  flow  below  the  vortex  sheet.  Again,  the  pressure  is 
lowest  at  the  erests  and  highest  at  the  troughs.  Since  the  crests  and  troughs  interchange  on  the 
two  sides  of  the  vortex  sheet,  the  result  is  that  a  net  pressure  imbalance  would  exist  across  the 
thin  mixing  layer  in  the  quasi-steady  approximation.  This  pressure  imbalance  is  in  phase  with  the 
vortex-sheet  displacement  and  hence  would  tend  to  increase  its  amplitude  leading  to  the  well- 
known  Kelvin-Helmholtz  instability. 


c  >  speed  of  sound 


Figure  3.2.  Pressure  distribution  on  the  two  sides  of  a  vortex  sheet  layer  at 
supersonic  convective  Mach  numbers  as  viewed  in  the  wave  frame  of  reference. 
+,  high  pressure  region;  -,  low  pressure  region. 


Now  suppose  the  Mach  number  of  the  flow  is  highly  supersonic  so  that  in  the  wave  frame  of 
reference  the  flows  on  the  two  sides  of  the  vortex  sheet  are  supersonic  as  shown  in  Figure  3.2. 
For  supersonic  flow  over  a  wavy  wall  it  is  known  that  the  pressure  distribution  is  no  longer  in 
phase  but  rather  90°  out  of  phase  with  the  displacement  of  the  wall.  This  results  in  identical 
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pressure  distribution  on  both  sides  of  the  vortex  sheet.  The  net  effect  is  that  the  wave  becomes 
neutrally  stable.  This  is  in  agreement  with  the  prediction  of  Miles  (1958)  using  hydrodynamic 
stability  theory. 


(a) 


c  <  speed  of  sound 


Figure  3.3.  Pressure  distribution  on  the  outside  surface  of  a  cylindrical  wavy 
vortex  sheet  jet  as  viewed  in  the  wave  frame  of  reference.  Also  shown  is  the  Mach 
wave  system  inside  the  jet.  (a)  Subsonic  flow  outside,  (b)  Supersonic  flow 
outside. 

For  circular  jets  the  situation  is  somewhat  different.  The  cylindrical  vortex  sheet  that  bounds 
the  jet  also  tends  to  reflect  acoustic  disturbances  that  impinge  on  it.  Thus,  acoustic  disturbances 
could  be  trapped  inside  the  jet  bouncing  back  and  forth,  forming  a  periodic  Mach  wave  system  as 
shown  in  Figure  3.3.  The  condition  under  which  such  a  Mach  wave  system  can  exist  is  that  in  the 
wave  frame  of  reference,  the  flow  inside  the  jet  is  supersonic.  Since  acoustic  disturbances  can 
propagate  upstream  (the  phase  velocity  c  of  the  wave  is  negative  in  this  case)  or  downstream  (c 
is  positive)  relative  to  the  flow  of  the  jet,  this  type  of  Mach  wave  system  exists  in  supersonic  as 
well  as  subsonic  jets  as  long  as  uj-c  is  greater  than  the  speed  of  sound  where  uj  is  the  jet 
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velocity.  Of  course,  for  subsonic  jets  c  would  have  to  be  negative,  namely,  the  wave  is  an 
upstream  propagating  wave.  The  pressure  distribution  associated  with  a  Mach  wave  system 
inside  a  cylindrical  wavy  wall  can  easily  be  calculated.  It  is  easy  to  show  that,  depending  on  the 
wavelength,  the  pressure  distribution  may  be  180°  out  of  phase  with  the  radial  vortex-sheet 
displacement.  The  wave  speed  c  may,  however,  be  subsonic  or  supersonic  relative  to  the  ambient 
gas.  Let  us  first  consider  the  case  of  subsonic  waves.  In  the  wave  frame  of  reference,  the  ambient 
flow  is  subsonic.  The  pressure  distribution  associated  with  subsonic  flow  outside  a  wavy  cylinder 
is  180°  out  of  phase  with  the  radial  displacement  of  the  cylindrical  wavy  wall.  Thus,  by  suitable 
choice  of  the  wavelength  of  the  Mach  wave  system,  pressure  balance  on  the  two  sides  of  the  thin 
mixing  layer  of  the  jet  is  possible  (see  Figure  3.3(a)).  This  implies  that  the  vortex-sheet  jet  can 
support  a  family  of  neutral  waves.  For  supersonic  jets  it  can  be  shown  that  if  the  effect  of  finite 
mixing-layer  thickness  is  included,  this  family  of  waves  is  actually  unstable.  Clearly,  with 
subsonic  phase  velocity  relative  to  the  ambient  gas,  the  amplitudes  of  these  waves  must  decay 
exponentially  in  the  radial  direction  outside  the  jet.  In  other  words,  the  disturbances  associated 
with  this  family  of  waves  are  confined  mainly  inside  the  jet.  These  subsonic  waves  are  the  third 
family  of  waves  observed  by  Oertel  (1980).  To  distinguish  the  waves  of  this  family  from  the 
other  two,  there  will  be  referred  to  as  the  subsonic  (instability)  waves. 

For  very  high-speed  jets,  the  phase  velocity,  c,  of  the  Mach  wave  system  may  become 
supersonic  relative  to  the  ambient  gas.  In  this  case,  the  flow  is  supersonic  both  inside  and  outside 
the  jet  with  respect  to  the  wave  frame  as  shown  in  Figure  3.3(b).  The  pressure  distribution 
associated  with  a  steady  supersonic  flow  over  a  cylindrical  wavy  wall  can  be  easily  determined. 
On  comparing  the  pressure  distribution  formulae  inside  and  outside  the  cylindrical  vortex  sheet, 
it  is  evident  that  pressure  balance  is  impossible  regardless  of  the  choice  of  wavelength.  Hence, 
unlike  the  case  of  subsonic  waves,  no  neutral  waves  are  possible.  Further,  it  is  easy  to  show  that 
for  certain  wavelengths  a  pressure  imbalance,  which  is  in  phase  with  the  vortex-sheet, 
displacement  is  possible.  The  net  result  is  that  with  the  internal  Mach  wave  system  a  highly 
supersonic  jet  can  sustain  a  family  of  instability  waves  with  supersonic  phase  velocities.  With 
supersonic  phase  velocity  relative  to  the  ambient  gas,  the  instability  wave  will  generate  a  Mach- 
wave-like  near  field  as  discussed  by  Tam  (1971)  and  Chan  &  Westley  (1973). 

In  the  past,  numerous  studies  of  the  Kelvin-Helmholtz  instability  waves  of  compressible  jets 
have  been  carried  out.  References  to  some  of  these  works  can  be  found  in  Michalke  (1984),  Tam 
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&  Burton  (1984),  Zaninetti  (1986,  1987)  and  others.  Most  of  these  investigations,  however, 
focused  primarily  on  the  Kelvin-Helmholtz  instability  in  subsonic  and  low  supersonic  jets.  In  this 
primer,  a  vortex-sheet  jet  model  as  well  as  a  more  realistic  jet  model  with  continuous  velocity 
profile  and  finite  shear-layer  thickness  will  be  considered.  Experience  indicates  that  the  simpler 
vortex-sheet  model  can  usually  provide  reasonably  good  estimates  of  the  phase  velocity  of  an 
instability  wave.  But  for  the  purpose  of  calculating  accurately  the  growth  rate  of  the  wave,  a 
finite-thickness  jet  model  is  necessary. 

3.1.1  Vortex-Sheet  Model  of  High-Speed  Jets 


iiO.x.t) 


Figure  3.4.  Perturbed  motion  of  a  round  jet  bounded  by  a  vortex  sheet. 

Consider  a  supersonic  jet  of  velocity  uj  and  radius  Rj  bounded  by  a  vortex  sheet  as  shown  in 
Figure  3.4.  Fet  {r,d^)  be  a  cylindrical  coordinate  system  centered  at  the  axis  of  the  jet  with  thex- 
axis  pointing  in  the  direction  of  the  flow.  On  starting  from  the  linearized  continuity,  momentum 
and  energy  equations  of  a  compressible  inviscid  fluid,  it  is  straightforward  to  find  that  the 
pressure  associated  with  small-amplitude  disturbances  superimposed  on  the  mean  flow  inside 
and  outside  the  jet,  pi  and  po,  are  governed  by  the  wave  and  convective  wave  equation, 
respectively. 


ayp.  =  0 


(3.1) 


f  ^ 


(3.2) 
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where  Uo  and  aj  (subseripts  o  and  j  denote  physieal  quantities  outside  and  inside  the  jet)  are  the 
speeds  of  sound  outside  and  inside  the  jet.  Let  ^(6U,0  be  the  radial  displaeement  of  the  vortex 
sheet.  The  dynamic  and  kinematic  boundary  conditions  at  the  vortex  sheet  r  =  Rj  are 

Pi=Po^  (3-3) 


Po  ^  ’ 


(3.4) 


J_ 

~Pj 


J  dx) 


(3.5) 


It  is  straightforward  to  find  that  separable  solutions  of  the  above  equations  and  boundary 
conditions  that  also  satisfy  the  boundedness  condition  at  r  =  0  and  r  ^  oo  are 


i{kx-\-nO-(o() 


(3.6) 


^  2  2  1/2  2  2  2  1/2 

where  =  {k  -  o)  !  a^)  ,  ((a)  -  Ujk)  I  Qj-k)  .  The  branch  cuts  of  rjo  and  r]i  are  taken  to 

be 


-2^<arg;7„  <  2^,  0  <  arg  77.  <  ;z-. 


The  wave  number  and  angular  frequency  k  and  co  are  related  by  the  eigenvalue  relation 


D(a),k) 


( '  =  derivative). 


(3.7) 


statement  A:  Approved  for  public 
release;  distribution  is  unlimited. 


84 


AEDC-TR-06-3 


Equation  (3.7)  is  not  new  (see,  e.g.;  Tam,  1971;  Chan  &  Westley,  1973).  Solutions  of  (3.7)  or 
the  roots  of  D{a),k)  are  the  wave  modes  of  the  jet.  The  nature  of  these  wave  modes  will  be 
discussed  later. 


3.1.2  Finite  Thickness  Shear  Layer  Model 

In  a  real  jet  the  mean  velocity  and  density  are  continuous.  Experimentally  it  has  been  found 
that  the  flow  velocity  is  uniform  in  the  central  part  of  the  jet.  Surrounding  this  uniform  core  is  a 
mixing  layer  with  a  velocity  profile  that  can  be  closely  approximated  by  a  half-Gaussian  function 
(see,  e.g.;  Troutt  &  McEaughlin,  1982).  Thus,  the  mean  velocity  profile  in  the  core  region  of  the 
jet  will  be  taken  as 


'Uj  , 

Uj  exp 

1 

IN 

1  -Q 

_ _ _ 1 

r  <h  , 

9 

r>n. 


(3.8) 


in  the  present  calculation.  In  (3.8)  h  is  the  radius  of  the  uniform  core  and  b  is  the  half-width  of 
the  jet  mixing  layer.  The  parameters  h  and  b  are  related  by  the  condition  of  conservation  of 
momentum  flux 


J  /7M  r  dr  =  2  PjUjRj . 
0 


(3.9) 


The  mean  density  p  is  related  to  the  mean  velocity  u  by  the  Crocco’s  relation  (Prandtl  number 
is  assumed  to  be  unity). 

It  is  easy  to  show,  starting  from  the  linearized  equations  of  motion  for  an  inviscid,  non¬ 
conducting  compressible  fluid,  that  the  equation  governing  the  perturbation  pressure  p  is  (see, 
e.g.,  Tam  &  Burton,  1984) 

/{r, 6, X, t)  =  ^r) exp[z(kx:  +  n0-  (Ut)] ,  (3.10) 
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p  =  0, 


(3.11) 
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where  a  =  {yp^l~^'^ .  The  loeally  parallel  flow  approximation  has  been  invoked  in  deriving 
(3.11).  This  equation,  together  with  the  boundedness  eondition  at  r  =  0  and  r  ^  oo  form  an 
eigenvalue  problem  for  0=  ry  {k)  ox  k  =  k{a)).  The  eigenvalue  ean  be  determined  by  integrating 
this  equation  numerieally.  Details  of  the  numerieal  proeedure  may  be  found  in  Tam  &  Burton 
(1984)  and  will  not  be  elaborated  here.  To  initiate  the  iteration  eyele  of  the  numerieal  proeedure 
the  solution  of  the  vortex-sheet  jet  model  may  be  used. 

3.1.3  Existence  of  Three  Families  of  Wave  Solutions 

The  braneh  points  of  the  funetion  rfo  and  rji  of  (3.6)  for  a  given  real  value  of  k  in  the  eomplex 
(2>-plane  are, 


(o  =  ±ka^. 


CO  =  ka 


^  Uj  ^ 


(3.12) 


respeetively.  The  braneh  euts  of  these  funetions  are  shown  in  Figure  3.5.  The  signifieanee  of  the 
braneh  points  is  that  regions  of  the  eomplex  ry-plane  to  the  left  of  the  left-hand  braneh  point 
represent  waves  for  whieh  the  flow  has  supersonie  eonveetive  Maeh  number.  Conveetive  Maeh 
number  is  defined  as  the  Maeh  number  of  the  flow  measured  in  the  moving  frame  of  the  wave.  It 
is  straightforward  to  show  that  for  every  point  in  this  region  of  Figure  3.5(b)  the  supersonie 
eonveetive  Maeh  number  relationship  uj  -  cojk  >  Uj  holds  {co=  C0r  +  icof  where  subseripts  r  and  i 
denote  real  and  imaginary  parts.  Similarly,  for  every  point  to  the  left  of  the  left-hand  braneh 
point  of  irjo  in  Figure  3.5(a),  whieh  is  for  the  statie  environment  outside  the  jet,  the  inequality 
I  C0rlk\  >  Oo  applies. 

It  is  also  easy  to  show  that  points  in  the  region  of  the  eomplex  (u-plane  to  the  right  of  the 
right-hand  braneh  points  in  Figure  3.5  again  represent  waves  for  whieh  the  eonveetive  Maeh 
number  of  the  flow  is  supersonie.  In  this  ease  the  inequalities  are. 


k  -  Uj 


>aj. 


o’ 
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inside  and  outside  the  jet,  respeetively.  The  remaining  region  of  the  eomplex  (u-plane,  namely, 
the  vertieal  strip  between  the  two  braneh  points,  represents  waves  having  subsonie  conveetive 
Maeh  numbers. 


(a)  Im  (w) 


Figure  3.5.  Braneh  euts  of  (a)  i/^and  (b)  rji  in  the  eomplex  (u-plane  for  a 
supersonie  jet.  Shaded  areas  are  regions  with  supersonie  eonveetive  Maeh 
number. 


Aeeording  to  the  instability  meehanism  deseribed  in  the  previous  seetion,  supersonie 
instability  waves  exist  in  a  jet  only  if  the  eonveetive  Maeh  numbers  (for  the  flow  inside  and 
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outside  the  jet)  are  supersonic.  In  terms  of  the  complex  ru-plane,  this  is  possible  only  if  the  jet 
Mach  number  is  highly  supersonic  such  that  the  left-hand  branch  point  of  Figure  3.5(b)  lies  to 
the  right  of  the  right-hand  branch  point  of  Figure  3.5(a).  In  other  words, 

Uj>aj  +  ag.  (3.13) 

The  vertical  strip  between  the  two  branch  points  satisfies  the  supersonic  convective  Mach 
numbers  criterion. 

Now  for  a  vortex-sheet  supersonic  jet  the  instability  wave  modes  are  given  by  the  zeros  of 
the  dispersion  function,  D(a>,k),  of  (3.7).  To  locate  the  zeros  in  the  complex  (u-plane  for  a  given 
value  of  k,  the  following  grid-search  method  has  been  found  useful.  To  implement  this  method 
the  region  of  interest  in  the  (U-plane  is  first  subdivided  into  small  subregions  by  a  rectangular 
grid.  The  values  of  the  dispersion  function  D{a),k)  of  (3.7)  are  calculated  at  each  grid  point.  A 
plotting  subroutine  is  then  called  that  performs  a  two-dimensional  interpolation  of  this  set  of 
values  and  constructs  the  two  families  of  curves  Re(Z))  =  0  and  Im(Z))  =  0.  The  intersection  of 
these  curves  provides  a  first  estimate  of  the  locations  of  the  zeros  of  D.  These  values  are  then 
refined  by  applying  Newton’s  iteration  method. 

Figure  3.6  shows  a  typical  example  of  the  zeros  of  the  dispersion  function  D{ct>,k)  found  by 
the  grid-search  method  for  a  high  supersonic  jet.  In  this  example,  the  axisymmetric  waves  (n  =  0) 
of  a  Mach  number  4.0  cold  jet  are  considered.  The  wave  number  kRj  has  been  set  to  be  equal  to 
15.0  in  the  calculation.  It  is  evident  in  this  figure  that  there  are  three  families  of  zeros,  or  wave 
modes.  The  isolated  zero  with  the  largest  temporal  growth  rate  is  the  familiar  Kelvin-Helmholtz 
instability.  Immediately  below  the  Kelvin-Helmholtz  zero  is  a  family  of  zeros  lying  in  the 
vertical  strip  (with  Im(ry)  >  0)  between  the  right-hand  branch  point  of  irjo  and  the  left-hand 
branch  point  of  rji.  These  are  the  supersonic  instability  waves.  To  the  left  of  the  supersonic 
instability  waves  is  another  family  of  zeros.  These  zeros  lie  on  the  real  a>  -axis  between  the  two 
branch  points  okk  =  ±ao-  These  waves,  therefore,  have  subsonic  phase  velocity  relative  to  the 
ambient  speed  of  sound.  They  are  the  subsonic  waves.  For  non-axisymmetric  wave  modes  with 
n==  1,2,3,...,  maps  similar  to  Figure  3.6  have  been  constructed  by  the  grid-search  method. 
Again  they  exhibit  three  sets  of  wave  modes  resembling  those  in  Figure  3.6.  The  above  results 
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apply  to  hot  jets  as  well.  For  hot  jets,  the  supersonie  instability  wave  modes  exist  even  at  a  mueh 
lower  supersonie  Maeh  number.  When  eondition  (3.13)  is  not  satisfied,  only  the  Kelvin- 
Helmholtz  and  the  subsonic  waves  could  be  found.  Experimentally,  Oertel  (1979,  1980,  1982) 
found  that  his  second  set  of  waves  (the  W  waves)  exist  only  when  uj  >  aj  +  Ug,  independent  of 
Mach  number  and  the  type  of  gas  used  to  form  the  jet.  This  is  in  total  agreement  with  condition 
(3.13). 


Figure  3.6.  Zeros  of  the  dispersion  function  in  the  complex  ru-plane. 

Mj  =  4.0,  cold  jet,  kRj  =  15.0,  n  =  0  mode. - ,  Re  (D)  =  0; - ,  Im(£))  =  0. 

T,  right-hand  branch  point  of  irjo,  ft,  right-hand  branch  point  of  rji.  A,  Kelvin- 
Helmholtz  instability;  #,  Supersonic  instabilities;  ■,  Subsonic  wave  modes. 

3,2  Resonances 


Resonances  occur  when  there  is  a  build-up  of  intense  oscillations  at  a  discrete  frequency. 
This  is  possible  when  two  criteria  are  met.  First,  there  is  a  mechanism  to  trap  acoustic  energy  in 
the  system  at  the  particular  frequency.  Second,  there  must  be  a  source  of  energy  that  feeds  into 
the  system.  For  an  engine  test  cell,  the  normal  modes  of  the  system  allow  energy  to  be  kept  in  the 
system.  Energy  at  frequencies  other  than  the  normal  mode  frequencies  are  leaked  out  or  radiated 
away.  Thus,  the  normal  mode  frequencies  are  the  possible  resonance  frequencies.  But  in  order 
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for  resonance  to  occur,  energy  at  the  resonance  frequency  must  be  supplied  to  the  system.  In  the 
case  of  engine  test  cells,  a  good  source  of  energy  is  the  instabilities  of  the  jet  flows.  However,  jet 
instabilities  are  confined  to  restricted  bands  of  frequencies.  Hence,  engine  test  cell  resonance 
tones  can  be  observed  when  there  is  an  overlap  between  a  frequency  band  of  the  jet  instability 
wave  and  the  normal  mode  frequencies  of  the  test  cell.  Since  a  test  cell  may  have  infinitely  many 
normal  modes,  it  follows  that  resonance  does  not  occur  at  all  normal  mode  frequencies;  only  at 
those  that  lie  within  the  unstable  frequency  range  of  the  jet  flow. 

For  long  ducts,  the  dispersion  relation  of  the  duct  modes  is  given  by  (2.10)  for  a  circular  duct 
and  (2.11)  for  rectangular  ducts.  These  duct  modes  have  a  propagation  speed  equal  to  the  group 
velocity  dcoidk.  The  group  velocity  can  be  found  by  differentiating  the  dispersion  relation  with 
respect  to  k.  Thus,  for  a  circular  duct,  the  dispersion  relation  and  group  velocity  for  the  duct 
modes  are. 


+  2A/n  —  k  + 
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=  0, 


(3.14) 
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(3.15) 


The  duct  modes  with  zero  group  velocity  have  the  significance  that  they  will  not  propagate 
away  once  they  are  generated.  Energy  fed  into  these  modes  will,  therefore,  stay  locally.  Their 
amplitude  could  build  up  into  violent  oscillations  if  there  is  a  continuous  input  of  energy.  Thus, 
these  wave  modes  may  also  lead  to  resonance. 

By  setting  dcoidk  =  0,  it  is  easy  to  find  from  (3.15), 


k  = 


MqCO 

<i.(>-^«)' 


(3.16) 


Statement  A:  Approved  for  public 
release;  distribution  is  unlimited. 


90 


AEDC-TR-06-3 


However,  k  and  ry  of  a  duct  mode  are  related  by  the  dispersion  relation  (3.14).  Hence,  on 
substituting  (3.16)  into  (3.14),  it  is  simple  to  find  the  frequencies  of  the  zero  group  velocity 
resonances.  They  are. 


CO  = 


(3.17) 


Similarly,  for  long  rectangular  ducts,  the  zero  group  velocity  resonance  frequencies  are. 


2 


n  m 


2 


(3.18) 


3,3  Receptivity 


Figure  3.7.  Excitation  of  shear  layer  instability  waves  by  sound. 
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In  the  literature,  there  are  a  number  of  experimental  and  theoretieal  investigations  (e.g.; 
Ahuja  et  al.,  1982,  1984)  aiming  speoifieally  at  obtaining  a  fundamental  understanding  of  the 
proeesses  by  whieh  external  sound  waves  exeite  the  intrinsie  instability  waves  in  shear  flows. 
These  studies  were  motivated  by  possible  applieations  of  the  aequired  knowledge  to  separated 
flow  eontrol  by  sound  (Ahuja  et  al.,  1982,  1984)  tone  exeited  jets  and  shear  layers  (Tam,  1978; 
Ahuja  et  al.,  1982;  Beehert,  1983;  Lepieovisky  et  al.,  1982;  and  Tam  et  al.,  1985)  and  sound- 
instigated  boundary  layer  instability  and  transition  (Shapiro,  1977;  Murdoek,  1980).  The  purpose 
of  the  diseussion  here  is  to  provide  a  simple  physieal  pieture  of  the  exeitation  proeesses  or 
reeeptivity.  The  simple  physieal  pieture  deseribed  here  should  eomplement  the  mathematieal 
theory  developed  in  Tam  (1978),  Tam  and  Morris  (1985)  and  Tam  (1981). 

Let  us  eonsider  the  exeitation  of  the  instability  waves  of  a  two-dimensional  free  shear  layer 
by  sound  as  shown  in  Figure  3.7.  For  eonvenienee,  the  two-dimensional  free  shear  layer  may  be 
regarded  to  form  downstream  of  the  trailing  edge  of  a  flat  plate.  Under  normal  eireumstanees,  the 
nature  and  physieal  eharaeteristies  of  the  aeoustic  waves  and  the  shear  layer  instability  waves  are 
so  different  that  the  two  types  of  waves  do  not  internet  with  eaeh  other  at  all.  In  other  words,  in 
the  ease  of  Figure  3.7  one  expeets  that  part  of  the  ineident  aeoustie  waves  will  be  transmitted 
through  the  two-dimensional  free  shear  layer  and  part  of  it  will  be  refleeted  baek,  but  no 
exeitation  of  the  intrinsie  instability  waves  would  take  plaee.  In  order  that  eoupling  of  the  two 
waves  would  oeeur  as  the  sound  waves  pass  through  the  shear  layer,  eertain  physieal 
prerequisites  must  be  satisfied.  The  neeessary  eonditions  are  as  follows:  (1)  matehing  of  the 
frequeneies  of  the  ineident  sound  wave  and  the  exeited  instability  wave;  (2)  matehing  of  the 
phase  veloeities  (in  the  flow  direetion)  of  the  two  waves.  The  first  eondition  is  easy  to 
understand  sinee  it  is  natural  to  expeet  the  exeited  instability  wave  to  have  the  same  frequeney  as 
the  ineident  sound  wave.  To  understand  the  seeond  eondition  it  is  best  to  eonsider  an  observer 
traveling  with  the  instability  wave.  If  the  phase  veloeities  of  the  waves  are  different  than  in  the 
moving  frame  of  referenee  the  observer  will  see  the  sound  wave  travel  past  by  him.  During  half  a 
period  the  sound  wave  amplitude  is  positive  and  the  other  half  period  is  negative.  The  net  result 
is  that  over  a  period  of  oseillation  any  interaetion  between  the  two  waves  will  be  eaneelled  out  to 
zero.  The  situation  is  totally  different,  however,  if  the  two  waves  travel  with  exaetly  the  same 
phase  speed  in  the  flow  direetion.  No  eaneellation  would  take  plaee  and  the  two  waves  would 
have  a  long  time  to  eommunieate  and  to  internet  with  eaeh  other.  It  is  to  be  noted  that  the  phase 
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velocity  of  a  wave  is  equal  to  the  frequency  divided  by  the  wave  number.  Therefore,  necessary 
condition  (2)  may  be  restated  in  terms  of  the  wave  numbers;  that  is,  to  excite  the  instability 
waves  of  a  flow  by  sound  there  must  be  a  matching  of  the  wave  numbers  of  the  two  kinds  of 
waves. 

Let  us  now  pursue  these  necessary  conditions  further  to  see  under  what  circumstances  they 
can  be  fulfdled.  In  general,  the  pressure,  pu  associated  with  an  instability  wave  of  frequency  co 
(see  Figure  3.7)  may  be  represented  by 


pix,y,t)=  ffyy-'-"',  (3.19) 

where  a,  is  the  wave  number  and  J{y)  is  the  amplitude  distribution  across  the  shear  layer. 
Similarly  the  pressure  field  of  the  incident  sound  wave,  ps,  may  be  represented  by 

pix,y,t)=B^y^^^-‘^‘\  (3.20) 

where  is  the  wave  number  (in  the  x-direction)  and  B  is  the  amplitude  of  the  sound  wave.  Here 
g(y)  is  included  to  allow  for  possible  variation  of  the  incident  sound  wave  amplitude  and  phase  in 
the  y-direction.  Under  normal  conditions,  and  are  very  different,  and  thus  the  phase 
velocities  of  the  waves  are  also  very  different,  and  there  should  be  very  little  interaction  between 
the  two  waves.  Now  if  the  amplitude  B  of  the  sound  wave  is  not  a  constant  but  varies  rapidly  in 
the  x-direction,  the  situation  changes  drastically.  Instead  of  (3.20)  the  expression  for  the  sound 
pressure,  ps,  becomes 


p,{x,y,t)=  B{xy'-"'^)=  (3.21) 


In  (3.21)  B{k)  is  the  Fourier  transform  of  the  spatially  varying  amplitude  B{x).  As  can  easily  be 
seen,  the  sound  wave  may  be  regarded  as  a  superposition  of  many  wave  components  with  the 
wave  number  {k  +  as).  The  wave  number  spectrum  is  no  longer  discrete  (see  Figure  3.8(a))  as  in 
the  case  of  constant  amplitude.  It  is  broadband  as  depicted  in  Figure  3.8(b).  With  rapid  enough 
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amplitude  variation  the  wave  number  speetrum  of  the  ineident  sound  wave  eould  be  very  broad; 
broad  enough  to  overlap  the  wave  number  of  the  instability  wave.  When  this  oeeurs,  as  is  shown 
in  Figure  3.8(b),  the  neeessary  eonditions  for  the  exeitation  of  the  shear  layer  instability  wave  are 
met  and  strong  eoupling  of  the  two  types  of  waves  is  expeeted. 


Figure  3.8.  Broadening  of  wave  number  speetrum  by  amplitude  variation,  (a) 
Constant  amplitude;  (b)  strong  amplitude  variation. 

Consider  now  the  ease  of  tone-exeited  jets  as  shown  in  Figure  3.9.  Without  loss  of  generality, 
we  will  assume  that  the  ineident  aeoustie  waves  are  from  an  external  souree  outside  the  jet.  The 
arguments  presented  below,  however,  would  apply  to  internal  aeoustie  disturbanees  coming  out 
of  the  nozzle  exit  as  well.  On  emerging  from  the  nozzle  exit,  the  jet  flow  consists  of  a  uniform 
core  surrounded  by  a  mixing  layer  with  inflexional  velocity  profile.  Such  a  velocity  profile 
supports  a  band  of  instability  waves.  From  the  standpoint  of  the  inflexional  instability  waves,  the 
amplitude  of  the  incident  sound  wave  rises  rapidly  near  the  nozzle  exit  to  an  approximately 
constant  value  slightly  downstream  as  shown  in  the  lower  portion  of  Figure  3.9.  Thus,  there  is  a 
very  abrupt  amplitude  variation  of  the  forcing  function.  Such  rapid  changes,  as  is  well  known, 
produce  an  extremely  broadband  wave  number  spectrum.  This  spectrum  easily  overlaps  that  of 
the  instability  waves.  Hence,  the  wave  number  matching  condition  is  fulfilled  and  strong 
excitation  of  the  instability  waves  by  the  incident  sound  waves  takes  place  in  the  immediate 
vicinity  downstream  of  the  nozzle  exit.  This  physical  reasoning  provides  readily  an  explanation 
to  the  well-known  experimental  fact  that  shear  layers  are  most  susceptible  to  acoustic  excitation 
downstream  of  a  trailing  edge. 
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Figure  3.9.  Forcing  function  as  seen  by  the  instability  waves  of  a  tone  excited  jet. 

In  the  above,  we  have  discussed  how  the  wave  number  matching  condition  can  be  satisfied 
when  there  is  a  rapid  amplitude  variation  of  the  incident  sound  waves.  This  is,  of  course,  not  the 
only  way  by  which  this  necessary  condition  can  be  satisfied.  It  is  easy  to  see  this  condition  can 
again  be  met  if,  instead  of  having  the  amplitude  of  the  incident  sound  wave  varying,  one  can 
have  a  rapid  spatial  change  in  the  mean  flow  velocity  profile.  In  flow  regions  where  this  is 
realized,  the  instability  wave  characteristics,  including  the  wave  number,  must  undergo  very 
rapid  spatial  changes  also.  Thus,  the  wave  number  spectrum  of  an  instability  wave  of  a  fixed 
frequency  is  broadband.  It  is  not  difficult  to  visualize  that  the  broadening  of  the  instability  wave 
number  spectrum  by  mean  flow  variation  could  be  sufficient  to  overlap  with  that  of  the  incident 
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sound  wave.  In  this  way  wave  number  matehing  ean  again  be  aeeomplished  and  strong  exeitation 
of  the  instability  wave  ensues.  This  type  of  flow  eondition  oeeurs  in  shear  layers  near  flow 
separation  regions  sueh  as  in  the  upper  surfaees  of  airfoils  at  high  angle  of  attaek.  The 
experiments  of  Abuja  et  al.  (1983,  1984)  are  exeellent  examples. 


3,4  Amplified  Instability  Wave  Spectrum 

A  high-speed  jet  is  unstable  only  over  a  finite  range  of  frequeneies.  As  the  instability  waves 
of  the  jet  are  the  energy  souree  of  resonanee  tones  it  follows  that  the  tone  frequeneies  must  he 
within  the  unstable  frequeney  band.  Linear  instability  wave  theory  indieates  that  the  growth  rate 
of  a  wave  depends  on  its  frequeney,  the  shear  layer  thiekness  or  the  mean  veloeity  profile,  the  jet 
Maeh  number  and  temperature,  and  the  azimuthal  mode  number.  Thus,  waves  at  different 
frequeneies  will  be  amplified  differently.  The  wave  with  the  largest  amplifieation  or  total  growth 
would  attain  the  highest  amplitude,  henee  it  would  most  likely  generate  a  tone  of  the  highest 
intensity.  Here  the  eharaeteristies  of  the  amplified  instability  wave  speetrum  are  diseussed. 

Linear  instability  wave  theory  is  now  a  relatively  well-established  subjeet.  The  mathematieal 
and  eomputational  proeedure  needed  to  determine  the  spatial  growth  rates  of  these  waves  are 
fully  doeumented  in  the  literature  (see,  e.g.,  Tam  et  al.,  1984,  1989). 

A  good  approximation  of  the  mean  flow  veloeity  profile  of  a  jet  is  given  by  (3.8)  and  (3.9). 
The  density  p  is  related  to  the  mean  veloeity  u  through  the  Croeeo’s  relation.  The  explieit 
relation  is. 


£lJ 

p  \ 


1  + 


Tr  K  T)u., 
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1 

2 

(3.22) 


where  Ta  and  7)  are  the  ambient  and  reservoir  temperatures,  respeetively. 

Insertion  of  the  mean  flow  quantities  into  (3.11)  leads  to  an  eigenvalue  problem.  The 
governing  equation  is  (3.1 1).  The  boundary  eonditions  are, 

r  =  0,  p  is  finite 
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r  ^  GO,  ^  ^  0  . 


For  a  given  frequency  instability  wave,  k  is  the  eigenvalue.  The  negative  value  of  the  imaginary 
part  of  k  or  -ki  is  the  local  growth  rate  of  the  instability  wave.  As  the  wave  propagates 
downstream  it  sees  a  continuous  change  in  the  mean  velocity  and  density  profiles.  Its  spatial 
growth  rate,  therefore,  changes  as  well.  For  a  wave  with  a  fixed  frequency / or  Strouhal  number, 
S  =fDjlUj,  and  azimuthal  mode  number  n  in  a  supersonic  jet  of  Mach  number  Mj  and  temperature 
ratio  TrITa,  the  total  amplification  is  given  by 


/ 


AyS,n,M^,W 


Tj 


(3.23) 


where  Xc  is  the  location  at  which  the  wave  reaches  its  maximum  amplitude  or  the  local  growth 
rate  is  zero.  Extensive  experimental  data  indicate  that  b  varies  nearly  linearly  with  x.  That  is, 
dbldx  =  cr  is  a  constant.  By  means  of  the  linear  relationship  between  b  and  x,  the  total  growth 
integral  on  the  right  side  of  (3.23)  may  be  rewritten  in  the  dimensionless  form  (except  for  a 
factor  a  ) 


1  = 


h 


R 


(3.24) 


where  bo  is  the  half  width  of  the  mixing  layers  at  the  nozzle  exit  and  be  =  h(Xc). 


statement  A:  Approved  for  public 
release;  distribution  is  unlimited. 


97 


AEDC-TR-06-3 


Figure  3.10.  Total  growth  integral  of  the  instability  waves  of  a  Mach  1.1  (cold) 

supersonic  jet  as  a  function  of  Strouhal  number  {/Dj/Uj): - n  =  0;  - 

n=  and ■  -n  =  2  mode. 


I 


Figure  3.11.  Total  growth  integral  of  the  instability  waves  of  a  Mach  1.4  (cold) 

supersonic  jet  as  a  function  of  Strouhal  number: - n  =  0; - n=  and  - 

■  -  ■  -  n  =  2  mode. 
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Figure  3.12.  Range  of  the  most  amplified  axisymmetric  and  helical  instability 
wave  modes  of  cold  supersonic  jets 


Figure  3.13.  Total  growth  integral  of  the  instability  waves  of  a  Mach  1.3  jet  at 

TjITa  =  2.0  vs  Strouhal  number: - n  =  0;  -  n  =  1;  and  -  •  -  •  -  n  =  2 

mode. 


statement  A:  Approved  for  public 
release;  distribution  is  unlimited. 


99 


AEDC-TR-06-3 


Figure  3.10  shows  the  values  of  the  total  growth  integral  of  the  lowest  three  azimuthal 
modes,  n  =  0,1,2,  as  funetions  of  Strouhal  number  for  a  Maeh  1.1  eold  jet  (7j.=  Ta).  This  figure 
shows  that  the  axisymmetrie  mode  (n  =  0)  has  a  slightly  higher  total  growth  than  the  other 
modes.  The  maximum  amplifieation  oeeurs  at  Strouhal  number  0.5.  The  helieal  mode  {n  =  1)  has 
a  slightly  smaller  total  amplifieation  that  peaks  at  a  Strouhal  number  of  0.4.  It  is  found,  however, 
that  the  total  amplilioation  of  the  axisymmetrie  mode  beeomes  smaller  relative  to  the  helieal 
mode  as  the  jet  Maeh  number  inereases.  At  jet  Maeh  number  1.4  the  helieal  mode  has  the  highest 
amplifieation  of  all  the  modes,  as  ean  easily  be  seen  in  Figure  3.11.  This  means  that  at  low 
supersonie  Maeh  numbers  the  axisymmetrie  mode  dominates.  The  dominanee  gradually 
diminishes  as  the  Maeh  number  inereases.  At  about  Mj=  1.3  the  helieal  mode  beeomes  the 
dominant  mode.  Figure  3.12  shows  the  most  amplified  instability  wave  Strouhal  number  for  the 
axisymmetrie al  and  helieal  mode  instability  waves  as  funetions  of  the  jet  Maeh  number  (Tj.  =  To). 
Obviously  the  axisymmetrie  mode  wave  peaks  at  a  higher  frequeney  and  is  an  important  wave 
mode  only  up  to  a  Maeh  number  of  1.3.  The  helieal  instability  wave  peaks  at  a  lower  frequeney. 
It  is  the  most  amplified  wave  for  Mj>  1.3.  The  shaded  regions  in  Figure  3.12  represent  the 
parameter  spaee  in  whieh  the  unstable  waves  have  signifieant  total  growth  and  henee  are  most 
likely  to  be  observed. 

The  growth  rates  of  the  Kelvin-Helmholtz  instability  waves  are  strongly  affeeted  by  the  jet 
temperature.  As  the  jet-to-ambient  temperature  ratio  inereases,  the  total  growth  of  the 
axisymmetrie  mode  instability  waves  beeomes  smaller  and  smaller  when  eompared  with  the 
helieal  mode  waves.  Thus,  they  beeome  less  and  less  important.  Figure  3.13  shows  the  total 
growth  integrals  of  the  n  =  0,  1,  and  2  modes  as  funetions  of  Strouhal  number  for  a  Maeh  1.3  hot 
jet  at  Tj/Ta  =  2.0  (7j  is  the  fully  expanded  jet  temperature).  Clearly,  the  helieal  instability  wave 
mode  (n  =  1)  is  dominant.  The  most  amplified  instability  wave  now  has  a  Strouhal  number  of 
0.13.  This  is  substantially  lower  than  that  of  a  eold  jet.  In  general,  the  Strouhal  number  of  the 
most  amplified  instability  wave  deereases  with  an  inerease  in  TjITa.  In  addition,  it  deereases  with 
an  inerease  in  jet  Maeh  number. 
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3.5  Super  Resonance 

As  previously  discussed,  a  duct  or  a  diffuser  can  support  a  variety  of  normal  duct  modes  of 
discrete  frequency.  If  one  of  the  frequencies  falls  into  the  unstable  range  of  a  jet,  then  resonance 
will  result.  Figure  3.14  shows  a  feedback  loop  formed  inside  a  duct  between  an  acoustic  normal 
mode  and  a  large  scale  instability  wave  of  the  jet.  The  jet  shear  layer  near  the  nozzle  exit  is  most 
receptive  to  excitations  induced  by  the  normal  mode  of  the  duct.  Since  the  shear  layer  of  the  jet 
is  thin  just  downstream  of  the  nozzle,  the  growth  rate  of  the  instability  wave  is  large.  Thus,  the 
excited  instability  wave  grows  rapidly  as  it  propagates  downstream.  When  the  instability  wave 
acquires  a  large  amplitude,  its  oscillatory  motion  feeds  energy  into  the  duct  mode.  The  duct 
mode,  in  turn,  excites  more  instability  waves.  In  this  way,  the  feedback  loop  is  closed  and  an 
acoustic  resonance  occurs. 

Now  as  pointed  out  in  the  last  section,  different  frequencies  of  an  instability  wave  mode 
attain  different  total  growth.  There  is,  however,  a  frequency  for  which  the  total  growth  is  the 
largest.  Now  if  the  frequency  of  the  duct  mode  in  the  feedback  loop  coincides  or  very  nearly 
coincides  with  the  frequency  of  maximum  growth  of  the  same  mode  of  instability  wave,  then  the 
feedback  would  be  most  intense.  Such  a  resonance  is  referred  to  as  super  resonance.  Super 
resonance  is  violent  and  could  lead  to  structural  damage  to  an  engine  test  cell. 


J.  J  /  /  .1  /  /  //  //  //  //  //  //  //  //  //■/. 


Figure  3.14.  Schematic  diagram  of  a  feedback  loop  formed  inside  a  duct  between 
an  acoustic  normal  duct  mode  and  a  large-scale  instability  wave  of  the  jet. 
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Part  II.  Recommendations 

1.  Suppression  of  Super  Resonance:  the  Method  of  Barriers 

Super  resonance  is  extremely  violent.  Once  it  occurs,  it  must  be  suppressed;  otherwise  the 
test  cell  could  suffer  severe  structural  damage.  Engine  development  may  be  interrupted.  In  the 
past,  AEDC  had  found  that  the  spraying  of  large  quantities  of  water  could  damp  out  the 
resonance  to  allow  the  testing  to  continue.  The  installation  of  Helmholtz  resonators  tuned  to  the 
resonance  frequency  also  was  somewhat  effective  in  reducing  resonance  amplitude. 

Here  as  a  part  of  the  consultant’s  report,  a  possible  new  method  for  the  suppression  of  super 
resonance  is  described  below. 

Super  resonance  is  driven  by  a  feedback  loop.  So  a  natural  way  to  suppress  such  a  resonance 
is  to  break  up  the  feedback.  The  feedback  loop  consists  of  two  major  components.  They  are  the 
instability  waves  of  the  jet  flow  and  the  normal  acoustic  modes  of  the  diffuser/duct  of  the  engine 
test  cell.  Now  for  high  temperature  circular  jets,  the  most  amplified  instability  waves  have 
azimuthal  wave  number  equal  to  ±1.  These  are  helical  or  flapping  mode  instability  waves.  Eor 
rectangular  jets,  the  most  amplified  waves  belong  to  the  antisymmetric  flapping  mode.  Now,  in 
order  for  the  instability  wave  to  couple  to  a  duct  acoustic  mode,  the  duct  mode  must  also  be  a 
flapping  mode  (or  a  combination  of  left  and  right  hand  helical  mode,  n  =  ±1).  This  is  a  necessary 
condition  for  coupling  the  oscillatory  motion  of  the  duct  mode  to  the  flapping  motion  of  the 
instability  wave  of  the  jet.  So  a  possible  suppression  method  is  to  insert  barriers  inside  a 
diffuser/duct  to  disallow  any  flapping  duct  modes.  This  can  be  achieved  by  installing  barriers 
inside  the  duct  so  that  there  is  no  symmetric  plane  across  the  center  of  the  duct.  A  flapping  duct 
mode  requires  such  a  symmetry.  A  barrier  configuration  that  disallow  flapping  duct  modes  is  as 
shown  in  Eigure  A.  It  is  easy  to  see  for  a  barrier  configuration  as  shown,  there  is  no  symmetric 
plane  through  the  center  of  the  duct.  In  the  absence  of  mode  matching  between  the  most 
amplified  instability  waves  and  duct  modes,  the  feedback  loop  is  broken  and  super  resonance 
may  be  avoided. 
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Figure  A.  Schematic  diagram  of  barriers  inserted  inside  a  diffuser  to  disallow  flapping  duct 

modes 

II.  Establishing  a  CAA  Capability  at  AEDC 

In  planning  for  the  establishment  of  a  CAA  capability  at  AEDC,  it  is  important  to  recognize 
that  this  capability  has  to  be  consistent  with  the  overall  mission  of  AEDC.  With  this  as  the 
background,  it  is  easy  to  identify  the  two  main  purposes  for  establishing  a  CAA  capability  at 
AEDC.  They  are; 

1.  To  support  engine  testing  and  development. 

2.  To  support  future  engine  test  cell  design  and  operation. 

The  principal  activities  associated  with  the  first  purpose  are: 

a.  Computation  of  the  acoustic  characteristics  of  existing  engine  test  cells. 

b.  Computation  of  the  aeroacoustic  characteristics  of  jet  engine  exhaust. 

c.  Offering  a  trouble-shooting  capability  for  unexpected  problems  encountered  in  routine  engine 
tests. 
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To  accomplish  (a)  and  (b)  engineers  may  use  eomputer  eodes  developed  in-house  or 
externally.  The  use  of  external  eodes,  however,  still  requires  a  eapable  in-house  engineer  to  run 
the  eode.  Short  term  problems  would  best  be  handled  by  an  in-house  CAA  person.  Longer  term 
problems  may  be  resolved  with  assistanee  from  outside  eonsultants. 

To  support  future  engine  test  eell  design  and  operations  would  require  an  in-house  engineer 
with  deeper  understanding  of  aeroaeousties  of  engine  test  eehs  and  CAA  methodology.  This  is 
neeessary  and  desirable  as  test  eell  geometry  and  operating  eonditions  may  ehange  during  the 
design  phase.  CAA  eodes  required  may  be  developed  in-house  with  eonsultant  assistanee  or 
developed  externally  with  input  and  partieipation  of  AEDC  engineers. 

PERSONNEL  REQUIREMENTS  —  It  is  believed  that  AEDC  should  have,  at  least,  one  on-site 
engineer  with  extensive  aeoustie  experienee  and  knowledge  of  CAA  methodology.  It  is 
important  to  point  out  that  just  having  programming  skill,  even  with  CED  experienee,  is  not 
sufheient.  Many  aeroaeoustic  problems  eannot  be  solved  by  eomputing  alone.  There  is  a  need  to 
know  or  have  a  fair  idea  of  the  physies  of  the  problem.  It  would  be  good  to  add  a  junior  CAA 
person  so  as  to  form  a  basie  team.  This  may  be  a  CED  person  with  strong  desire  to  switeh  over  to 
the  CAA  area.  They  ean  work  with  existing  CED  engineers  to  provide  a  eore  CAA  eapability  at 
AEDC. 

CODE  DEVEEOPMENT  STRATEGY  —  In  CED,  an  often  adopted  strategy  is  to  aequire  an  all¬ 
purpose  CED  eode.  Sueh  a  eode  would  allow  engineers  to  do  numerous  types  of  CED  problems. 
This  is  a  good  strategy  for  CED,  as  most  engineers  would  only  need  to  learn  to  run  a  single  eode. 
In  CAA,  this  may  not  be  the  best  strategy.  The  reason  is  simple.  CAA  problems  are  often  very 
dissimilar.  Eaeh  elass  of  problems  has  its  own  oharaeteristies.  As  a  result,  there  is,  at  the  present 
time,  no  all-purpose  CAA  eode.  Also  the  run  time  for  CAA  eodes  is,  in  most  oases,  very  long.  To 
keep  the  run  time  reasonable,  it  is  neeessary  to  use  eodes  that  are  optimized  for  the  partioular 
elass  of  problem  one  wishes  to  solve.  All-purpose  eode  eannot  be  highly  optimized.  So  the  run 
time  will  inevitably  be  fairly  long.  With  the  above  in  mind,  it  is  reoommended  that  the 
development  of  a  set  of  dedioated  CAA  eodes  for  engine  testing  and  development  purposes  be 
regarded  as  a  priority  in  the  establishment  of  a  CAA  eapability  at  AEDC. 
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